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INTRODUCTION 


The surface of the earth continually atforbs and emits energy. 

The balance achieved between energy inputs and energy outputs represents 
the energy budget. Models have been developed to simulate the complex 
interactions among the energy budget components. Ground-based meteoro- 
logical point-sample measurements have been used as inputs for verifi- 
cation of model accuracies. These models can be used with data from 
aircraft and satellites equipped with sensors capable of measuring 
reflected and emitted energy. The remote sensor offers a synoptic, 
yet detailed, approach for studying and using concepts from the energy 
budget for resource monitoring over a variety of landscape variables 
on a repetitive basis. 

The temperature of a fallow soil or vegetated surface is an integral 
response by that surface to environmental factors such as level of 
incoming radiation, air temperature and humidity, wind speed, etc. 

Such surface characteristics as albedo, roughness, wetness, and other 
physical parameters alter the effect of the environmental factors on 
the energy budget and heat exchanges and consequently are expressed 
as temperature variations. 

A major source or sink of heat energy in the soil -pi ant environment 
is water. Water plays a major role in conduction and convection of 
heat to and from the soil surface. In addition, phase changes of water 
by vaporization and fusion absorb and release large quantities of heat. 
Fallow soil surfaces dissipate a portion of the incoming energy by 
evaporation. Dissipation of energy by evaporation from plants is termed 
transpiration. A land surface having a crop canopy dissipates energy 
by evaporation from both the plant and the soil which is collectively 
termed evapotranspi ration (ET). Soil moisture and evapotranspiration 
estimates using remote-sensor data could provide a powerful management 
tool for irrigation scheduling, crop yield predictions, detecting 
disease and insect infestations, flood forecasting, predicting seed 
germination, monitoring land-erosion potential, and numerous other 
appl i cations. 
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OBJECTIVE 

To determine the physical and thermal properties of the surface 
and subsurface soils and to measure the surface temperature and emitted 
radiation by use of remotely sensed data. 

BACKGROUND 

Energy Budget 

According to Sellers (1965), approximately 30% of the solar radiation 
reaching the outer boundary of the earth's atmosphere is reflected 
and scattered back into space by clouds, aerosols, and other atmospheric 
constituents. Approximately 17% of the incoming radiation is absorbed 
by the earth's atmosphere and about 22% reaches the earth's surface 
as diffuse radiation. Therefore, only about 31% of the radiation entering 
the outer boundary of the earth's atmosphere reaches the surface as 
direct solar radiation. 

Incoming radiation received at the earth's surface is either ab- 
sorbed, reflected, or transmitted. Absorbed short-wave radiation may 
be emitted at a later time at a longer wavelength; therefore, solar 
radiation may be absorbed, stored in the crop canopy or soil, and emitted 
later as thermal radiation. A portion of the absorbed radiation can 
be used to evaporate water. The heat of vaporization used in the evapo- 
ration process may subsequently be transferred from the surface as 
a vapor flux. 

That portion of the radiation in the energy budget which is the 
difference between total incoming and total outgoing radiation is termed 
net radiation. Mathematically, net radiation (R^) can be expressed as 

'*n " (‘’s * '’s 

where R^ is direct incoming shortwave radiation, R^ is diffuse incoming 
shortwave radiation, Rj is incoming longwave radiation, Rfo Ts outgoing 
longwave radiation, and a is the albedo of the surface. Since albedo 
directly affects the net radiation, remote monitoring of albedo provides 
information about the energy absorption characteristics of a surface. 

The absorbed energy is potentially available for the evaporation of 
water, for heat storage, etc. In addition, reflectance within specific 
wavebands may serve as an indicator of certain physical and biological 
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properties. For example, surface water has a unique spectral signature 
in the longer reflective- infrared wavelengths which can be used to 
thematically map the occurrence of water. The evaporation rate of 
free water surfaces is limited only by energy budget conditions and 
not by the transport of water to the evaporating surface. 

The net radiation may be utilized at the earth's surface as 

R+S+A+LE+P+M=0 [2] 

where S is soil heat flux, A is sensible heat flux in the air, LE is heat 
flux due to evapotranspiration, P is the energy used in photosynthesis, 

and M represents a miscellaneous energy term. For some applications, 
equation [2] is used as an energy balance equation where only surface 
elements are considered and storage plus divergence terms are neglected. 
However, energy may be transferred into or away from a volume element 
of the surface by advection of sensible and latent heat which may be 
expressed as 

Advection = Cp v(pyT) dz + ^ dz [3] 

where z is distance along the path, ? is a defined upper boundary, C 

3 3 ^ 

is the specific heat of air, v is -^ + p is air density, y is wind 

speed, T is temperature, L is latent heat of vaporization, e is water 

vapor flux, R is the universal gas constant, and e is the vapor pressure 

of the air. 

Energy can be stored in a volume element of the earth's surface 
in the form of heat content of the crop, sensible heat in the air, 
and latent heat of the air. These storage terms may be represented as 

storage = Qpc (|I) dz + C^p (|I) dz t ^ (||) dz [4] 

where Q is the heat capacity of the crop, pc is the crop density, e 
is the ratio of molecular weight of water to air, t is time, and the 
other terms are as previously defined. 

Utilizing the energy budget terms developed for an element 6x6ySz 
of the earth's surface when a crop canopy or other rough surface is 
present, an energy budget equation can be formulated. The complete 
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energy balance equation is 

+ S + A + LE + P + M + /^ Cp v(pyT) dz + ^ v{^) dz + 

Qpc (|I) dz + /' Cpp (|I) dz + ^ (||) dz = 0 [5] 

where all terms have been previously defined. 

The relative magnitude of the storage terms is negligible during 
most of the day. However, the error introduced by neglecting the 
storage terms can be appreciable for a short period near sunrise when 
A, and E are small and |^is large (Tanner, 1960). Tanner found 
that at night (midnight to 0600 hours) the storage terms of an alfalfa 
brome canopy were about 6% of the soil heat flux and about 1 % of the net 
radiation. The horizontal divergence terms may be large if the area is 
small and permits air passage through the crop (Tanner, 1960). Divergence 
may also be high where an irrigated field is surrounded by a dry area. 
Divergence is normally negligible well inside large fields where border 
or other external influences are small. 

The flux direction of the energy balance components are dependent 
on time of day and climatic conditions. During the day, radiation 
from the sun is absorbed by the soil and the crop. Thus, the soil 
and crop surfaces may become warmer than the air so that sensible heat 
flux is away from the crop volume. During the day under stable con- 
ditions, the soil surface is warmer than the soil at greater depths 
so that the soil heat flux is downward away from the crop volume. 

At night, soil and crop surfaces lose heat through the emission of 
longwave radiation which is nearly always greater than atmospheric 
counter-radiation; therefore, soil heat flux is toward the crop volume 
to counteract the radiation loss. At night, the air is normally cooler 
than the crop and soil surfaces resulting in sensible heat flux away 
from the crop volume. 

If water is available for evaporation, the latent heat flux is 
away from the crop volume during the day. At night, latent heat flux 
may be away or toward the crop volume depending on the respective 
temperatures of the air, crop, and soil surfaces. In regions with 
dry climates, latent heat flux will normally continue to be away from 
the crop volume both day and night and dew formation will not occur. 
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However, when the crop or soil temperature is less than the dew point 
temperature, the latent heat flux will be toward the crop canopy. 
Horizontal divergence terms will normally be greater during the day 
than during the night due to higher wind speed. 

The energy balance components which have been discussed are for 
micro-climatic regimes such as agricultural fields or portions of fields. 
Over large macro-climatic land areas, horizontal divergence terms may 
be neglected. For long-term periods, storage terms may be neglected 
because energy stored at one time is ultimately released at another 
time. For example, heat stored in the soil and crop during the day 
is released at night. 

Evapotranspi ration 

Equation [5] may be used to describe the energy balance of a 
vegetated land surface. If advection, storage, and miscellaneous energy 
terms are neglected, equation [5] becomes 

Rn + S + LE + A = 0 [6] 

for a cv’op canopy. The energy balance equation caij be solved to determine 
evapotranspi ration utilizing the Bowen ratio (Bowen, 1926). The Bowen 
ratio (B) is the ratio of sensible heat (A) transport to latent heat 
(LE) transport. Using the Bowen ratio, equation [6] may be written as 

LE = -(Rn + S) / 1 + B [7] 

where LE is the evaporative flux or evapotranspiration rate for a 
cropped surface and B is the Bowen ratio or A/LE, Additional predictive 
models using energy budget, mass transfer, or combination equations 
have been developed to estimate the evaporative flux. Bartholic, Namken, 
and Wiegand (1970), McGuinness and Bordne (1972), Stone and Horton 
(1974), and Rosenberg, Hart, and Brown (1968) have reviewed and evaluated 
several of the models. The prediction normally estimates the potential 
ET in contrast to the actual ET. Morton (1969) has stated that with 
adequate moisture supplied to the evaporating soil and vegetated surfaces, 
the actual ET is accurately estimated from the potential ET on a regional 
basis. However, this relationship is erroneous for surfaces which 
are not well watered. Remote sensing may provide a method of detecting 
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and mapping these well -watered and dry surfaces far improving the 
application of existing models. 

An example ET equation is the one derived by Wiegand and Bartholic 
(1970) which uses surface temperature in estimating potential evaporation 
from a wet surface. The crop canopy is assumed to be a wet surface 
at the saturation vapor pressure; therefore, the canopy vapor pressure 
is a function of the prevailing canopy temperature. Evapotranspiration 
as expressed by their equation is 

Ep - -(Rn+S) / {1 + G [(VT^)/e; - e^)]} [8] 

where Ep is the potential evapotranspiration, is air temperature 
at height 'a' above the surface, is surface temperature, e^ is the 
saturated water vapor pressure at height 'a', and e^ is the saturated 
vapor pressure of the evaporating surface at temperature T^. Equations 
employing a surface temperature parameter can be applied for evapotrans- 

4 

piration predictions using remotely sensed emittance estimates of tem- 
perature. However, when a surface is comparatively dry or the soil 
moisture is limiting to the crop canopy, the air in the canopy is not 
at saturation so the vapor pressure can not be predicted using the 
remote temperature measurements. 

The Jensen-Haise model for potential ET prediction (based on tem- 
perature and radiation) was developed and evaluated over large areas 
of the world (Jensen and Haise, 1963). The model is extensively used 
for irrigation scheduling and has input parameters as 

ET^ = (0.025 T + 0.08) R [9] 

p s ^ 

where ETp = potential ET in mm min"' 

T = air temperature in degrees C 

= equivalent depth of evaporation of incoming solar 
radiation in mm/mi n"^. 

Although air temperature is used, the temperature of the evaporating 
surface should be an improved parameter for estimating actual ET. 

The same is probably true for R^. 

The Jensen-Haise ETp is frequently modified for actual evapotrans- 
piration using a crop coefficient, Kc. The crop coefficient varies 
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from near zero to slightly more than 1.0, depending upon the growth 
stage and type of crop. During the period of maximum vegetative cover, 

Kc will be a maximum for that crop. An apparent crop coefficient may 
be derived from remotely sensed data using reflectance and emittance 
properties of the vegetated surface in contrast to the present tedious 
methods of obtaining crop coefficients. 

Water Budget 

A water budget is frequently used to estimate water used in evapo- 
transpiration over a specific time interval. The water conservation 
equation for a volume of plant root zone for a specific time period is 

P+I = R + aW + E + T + D [10] 

where P is the precipitation, I is irrigation, R is surface runoff, 

AW is the change in stored water in the soil volume during the time 
interval, E is evaporation from soil, T is transpiration from the plants, 
and D is the amount of soil water either entering or leaving the soil 
volume. The precipitation, irrigation, runoff, and change in storage 
terms are readily measurable. The profile drainage, evaporation, 
and transpiration terms are more difficult to measure. 

Between rainfall or irrigation applications, the terms P, I, and 
R will equal zero. The water conservation equation then becomes 

-AW = E + T + D [11] 

with the terms previously defined. The change in water storage (aW) 
is measured as profile water content at the end of the time period 
minus profile water content at the beginning of the period. Equations 
(10) and (11) are both in the integral form with the terms being totaled 
over a given period of time. It is often desired to speak of water 
loss as a rate, or loss per unit of time. Equation (11) expressed 
in differential form becomes 

-dW/dt = dE/dt + dT/dt + dD/dt [12] 

where W is the water storage in the soil volume and t is time. Equation 
(12) is the time rate of change form of the water conservation equation 
between water applications. The rate of evaporation from soil (dE/dt) 
and the rate of transpiration from plants (dT/dt) are combined to form 
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the evapotranspi ration (ET) rate term. The rate of profile water 
drainage (dD/dt) is commonly referred to as the soil water flux (v). 

Equation (12) takes the form 

-(dW/dt) = ET + V [13] 

where (dW/dt) is the rate of water storage change in the soil volume. 

The flux term in equation (13) can be estimated using Darcy's 
equation 

V = Ki [14] 

where v is the soil water flux, K is the unsaturated hydraulic conduc- 
tivity, and "i" is the hydraulic gradient or the driving force. This 
computation requires considerable ground information, which is difficult 
to derive and for large areas is economically infeasible. 

When the water budget is used, neglect of terms describing water 
loss other than evaporative flux may introduce considerable error. 

For example. Black, Gardner, and Tanner (1970) determined on a field 
water budget study that during a 60-day study period the total water 
loss from a field of snap beans was 35 cm with 17 cm lost to ET and 
18 cm lost due to drainage from the 150-cm profile. Goltz et al . (1971) 
found that for an onion crop, the drainage loss exceeded the ET loss. 
Therefore, the water budget approach is only applicable where a consider- 
able knowledge of the landscape physical properties and the groundwater 
status is known. If the crop canopy and fallow soil surfaces are an 
indicator of the surface and subsurface moisture at the time of measure- 
ment, the errors in the traditional water budget would not be present 
where the estimation of soil moisture occurs at a point in time by 
remote sensors which would directly determine the stored soil water. 

Remote-Sensing Applications to Soil Water-Evaporation Assessments 

When the supply of water to the evaporating surface (primarily through 
conduction through the soil matrix to the fallow soil surface or translocation 
through plant roots and tops to the leaves) is limited, the evaporation 
rate decreases with a subsequent increase in temperature of the evaporating 
surface. In early work. Tanner (1963) used energy balance relationships 
to show that if the transpiration rate of plants decreased and if the 
radiation balance and wind structure remained the same, the decrease 
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in latent heat exchange resulted in an increased plant temperature. 

The incoming energy to the surface was utilized in heating the surface 
with dissipation by sensible and soil heat flux in contrast to being 
partitioned into latent energy with subsequent dissipation as latent 
heat flux. As moisture stress increases, with consequent decreases 
in ET rate and increases in surface temperature, the relative turgidity 
of the plant decreases. Wiegand and Namken (1966) observed in cotton 
plants that a decrease in relative turgidity of plants from 83 to 59% 
resulted in a 3.6 C increase in leaf temperature. This magnitude of 
temperature variation is well within the resolution capabilities of 
modern remote-sensing systems. Moore et al . (1974) illustrated that 
the X/5 thermal detector on the 5192 of SKYLAB had sufficient thermal 
resolution to detect emittance variations which were associated with 
soil moisture related terrain features in vegetated landscapes. 

Werner and Schmer (1972) investigated the application of remote 
sensing techniques to the inventory of soil water in croplands and 
rangelands. In addition to data from several film-filter combinations 
within the reflective region, thermal scanning data were collected. 

They concluded that the blue region of the electromagnetic spectrum 
could be used for determining soil water conditions with little or 
no crop cover. The green, red, and near-infrared regions were all useful 
for evaluating soil water where crop or range canopies were present. 

Their evaluations were conducted within specific agricultural land 
uses. The reflectance differences among crop types are considerable 
and vary with phonological stage; therefore, the reflective spectral 
region has limited use for predicting moisture variation within diverse 
agricultural regions. Since water has a considerable affect on the 
energy budget, surface emittance variations may provide a better estimate 
of soil moisture less influenced by crop cover. Idso, Jackson, and 
Reginato (1975) have presented a comprehensive discussion of the detection 
of soil moisture on fallow land. They illustrated the use of a predicted 
moisture tension in contrast to absolute moisture quantity to reduce 
the variation among soil types. 
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TEST SITE IDENTIFICATION 
Test-site Selection Procedure 

Many test sites were selected for anticipated data collection 
during the SL-2, SL-3, and SL-4 missions. Site selections were made 
by the investigators and personnel familiar with the area by evaluation 
of historical LANDSAT and aircraft imagery and by field reconnaissance 
surveys. After initial site selection, a general field survey was 
conducted by on-site ground inspection. The criteria for site selection 
included the following: uniform soils, a diversity of agricultural 

land uses, the occurrence of irrigated crop and fallow fields with 
dryland of similar land use in adjacent areas, and landowner approval. 
Since energy budget information was of primary interest, thermocouples 
and heat flux plates were placed in the soil at various depths to 100 cm. 
The minimum time required for implanting the sensors prior to measurement 
was 24 hours or one diurnal cycle. Therefore, considerable travel 
to potential test sites with detailed inspection was required, especially 
during SL-3 when the workshop was not in a stable orbit. 

Description of Carrizo Springs Test Site 

The site near Carrizo Springs, Texas (28° 30'N, 99° 50'W), was 
selected in accordance with the estimated coverage of the S-192 scanner 
and with the general criteria previously listed for test site selection. 
The soils were predominantly of the Uvalde-Montell -Atco Association 
and are nearly level to gently undulating, deep, silty clay loams, 
clays, and clay loams that have moderately and very slowly permeable 
subsoils (Soil Conservation Service, 1972). An image of SKYLAB S190-B 
color-infrared photography (Fig. 1) includes the test site and surrounding 
area. The predominant land use was range with some fields of irrigated 
agriculture, including alfalfa, cabbage, onions, carrots, spinach, 
and pre-i rrigated fallow. The dryland agriculture included fallow, 
improved range, and unimproved range. Thirty separate fields identified 
in Fig. 2 for which ground data were acquired were typical of the regional 
land-use and soil moisture variations. Area "9" was a surface water 
body. A frost had occurred approximately one-month prior to the SKYLAB 
overpass; therefore, the dryland agriculture was essentially dormant. 




hotoqraohy ac 
The two regi 
s. 2a and 2b 
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Fiq. 2a - Area A 



Color-infrared photograph illustrating the two intensive 
test sites near Carrizo Snrings, Texas. Data are from RB-57F 
NASA mission 260, RC-8 "6-inch" focal length camera at 18,300-m 
AGL. Refer to Table 1 for descriotion of land use and soil 
moisture for fields as numbered. Area "A" is a seep and areas 
"B" are center pivot irrigation systems (original in color). 
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DATA COLLECTION 

Ground Data 

Ground-based energy budget data were acquired for specific fields 
within the test site. Measurements included thermocouple soil temper- 
atures at several depths, air temperature and vapor pressure at several 
heights, net radiation, filtered radiometry, and soil heat flux. A 
description of the sensors and their placement is in Table A-1 of 
Appendix A. Only those measurements required for the specific evaluations 
in this document are presented. Land-use and soil moisture data were 
acquired to characterize the thirty agricultural fields. Proportions 
of green vegetation, fallow, and dry organic debris were determined 
from slides and are reported as percent of total area. Slides were 
exposed using a hand-held 35-mm camera vertically 1-m above ground 
level (AGL) over the test site. Measurements were conducted by projecting 
slides on a grid containing 100 randomly placed dots. Each value was 
reported as the mean of four separate observations on each of three 
replicate slides for each field. Soil moisture samples were collected 
in triplicate with each sample composed of two subsamples. The percent 
gravimetric soil moisture was determined for the 0->2 cm, 2->10 cm, and 
10^30 cm depths by oven drying at 105 C. A composite value for the 
30-cm profile was computed by weighting and summing the averages for 
the depths represented. 

The locations of the two primary areas of ground data collection 
are identified on the S190-B data in Fig. 1. The intensive test fields 
are numbered on prints of aircraft data presented in Fig. 2o Table 
1 is a summary of the ground data with field numbers corresponding 
to those in Fig. 2. 

Aerial Data 

A multistage data collection mission was accomplished including 
NC-130B, RB-57F, and SKYLAB data. The various operable sensors of 
each platform along with delivered data products are described in Table 
A-2 (Appendix A). The mission numbers include: 258 for NC-130B, 260 

for RB-57F, and orbit 94 on ground track 20 at 20:12 GMT (13:12 local 
CST) on 28 January 1974 for SKYLAB. The S-192 scanner aboard the SKYLAB 
workshop used the X/5 thermal detector for this pass. 
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TABLE 1. LAND USE AND SOIL MOISTURE OF AGRICULTURAL 
FIELDS WITHIN TEST SITE 


Field 

No.-^ Crop % Green-^ 

% Gravimetric 

Soil Moisture 

0-^2cm 

2->-l 0cm 

1 0->30cm 

Composite 

1 

range 

1.2 

11.33 

16.24 

15.06 

15.13 

2 

onions 

28.0 

11.96 

19.14 

18.50 

18.23 

3 

alfal fa 

22.8 

14.74 

18.50 

12.53 

14.27 

4 

alfalfa 

86.8 

35.21 

25.38 

22.94 

24.41 

5 

alfalfa 

86.1 

23.36 

19.97 

18.08 

18.94 

6 

fallow 

0.0 

8.72 

20.23 

24.81 

22.52 

7 

fallow 

2.6 

32.86 

31.98 

28.17 

29.50 

8 

fallow 

0.0 

12.42 

21.50 

19.00 

19.23 

9 

water 






10 

fallow 

6.0 

27.05 

26.67 

24.99 

25.57 

11 

alfalfa 

89.2 

31 .90 

23.67 

23.17 

23.91 

12 

alfalfa 

10.3 

15.63 

21 .09 

13.82 

15.88 

13 

fallow 

0.0 

11.40 

20.98 

18.26 

18.53 

14 

cabbage 

30.5 

5.81 

12.67 

14.77 

13.61 

15 

onions 

11.6 

3.77 

11.11 

12.50 

11.55 

16 

range 

1.2 

5.47 

9.04 

8.18 

8.23 

17 

cabbage 

10.3 

10.93 

14.83 

15.95 

15.32 

18 

cabbage 

82.1 

9.23 

13.00 

10.73 

11.24 

19 

cabbage 

85.0 

16.74 

20.32 

22.24 

21 .36 

20 

onions 

10.6 

9.20 

18.76 

19.89 

18.88 

21 

onions 

1.2 

6.78 

13.55 

14.80 

13.93 

22 

fal low 

0.0 

2.64 

10.95 

10.94 

10.39 

23 

carrots 

71.2 

11.15 

12.90 

13.86 

13.42 

24 

cabbage 

37.6 

22.66 

22.71 

22.13 

22.32 

25 

cabbage 

55.8 

12.44 

21.74 

21.59 

21.02 

26 

fallow 

0.1 

9.53 

22.29 

23.00 

21.91 

27 

fal 1 ow 

0.0 

5.95 

16.82 

13.85 

14.12 

28 

carrots 

62.0 

9.04 

18.11 

20.93 

19.39 

29 

carrots 

22.6 

4.93 

13.97 

16.14 

14.81 

30 

fallow 

0.2 

4.33 

14.56 

15.30 

14.37 

31 

range 

0.0 

7.01 

11.97 

10.55 

10.69 

A 

seep 






B 

center pivot irrigation 





— Refer to Figs. 2a and 2b for field locations. 


2 / 

— Reported as percent of total area which was green vegetation when 
viewing the surface vertically. 
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PROCEDURES FOR DATA REDUCTION 


General 

Visual and statistical analyses were conducted on the data. 

Visual analyses of the scanner data were possible using screening film 
products from the scanners and using the DAS analysis system at NASA, 

JSC. Statistical analyses of digital scanner data were conducted on 
unrectified products from the NC-130B and SKYLAB scanners. Film products 
from the RS-7 scanner aboard the RB-57F aircraft and the S190-A film 
from SKYLAB were digitized using the scanning microdensitometer {^ignal 
Analysis and dissemination d^uipment - SADE) at the Remote Sensing 
Institutec Location and statistical analyses of the thirty individual 
fields were conducted by: applying an edge detection algorithm to 

locate radiometrical ly homogeneous areas within the digital matrix, 
computing means and standard errors of the radiance values within 
each area, and testing on a field-by-field basis the field means for 
each spectral region using a t-test at the 0.05 confidence interval. 

Edge Detection 

Since the original S192 scanner data were in a conical form, rec- 
ognition and accurate mapping of field boundaries in the digital data 
were difficult. Therefore, an algorithm was developed to locate edges 
between relatively homogeneous areas in a matrix of digital data. 

The assumptions in the algorithm included: 


1. An ideal edge between two homogeneous data clusters occurs 
at a step change in the data. 

2. The rate of change of data values (gradient vector magnitudes) 
within a homogeneous area is small. 


3„ The gradient magnitude near and at an edge is greater than 
within a homogeneous area. 

4c Th.e gradient vector direction at an edge is perpendicular 
to an edgOo 

5o An edge point has the maximum gradient value when compared 
to adjacent points in the gradient vector direction. 


A general flow chart of the algorithm is presented in Fig. 3. 

The magnitude jF(x,y)| and the direction (e) are computed for each point by 


e = arctan (If /If ) [15] |F(x,y)| = (|f) 


3F 

9X 


1/2 


[16] 




Fig. 3 - Flow diagram of edge detection procedure which was 
applied to matrices of digital data. 
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sF aF 

The orthogonal derivatives and are approximated at the midpoint 
f four data points with the two numerical differences across the diag- 
onals of the points. The gradient is quantized into four discrete 
directions illustrated in Fig. 3. For example, direction one occurs 
when the gradient direction has a numerical value in the range of 22.5° 
to -22.5° or 157.5° to 202.5°. 

After initial computation of the gradient vector direction and 
magnitude, the edge detection algorithm is applied. For each separate 
data point, the gradient magnitude is compared to a threshold value. 

If the value is less than the threshold, the data point is not a potential 
edge point and the algorithm goes to the next point. If the gradient 
value is greater or equal to the threshold, the point is retained for 
further analysis. The gradient vector direction is perpendicular to 
an edge. An edge is located by comparing adjacent gradient vector 
magnitudes in the gradient vector direction. The four possible gradient 
directions are illustrated with data points a, b, and c in Fig. 3. 

An edge is located at point "b" when the gradient magnitudes at both 
adjacent points "a" and "c" in the gradient vector direction are smaller 
than the gradient magnitude at point "b". If this condition is not 
met the algorithm returns and goes to the next point. Edge points 
are coded with the gradient magnitude at that point and other non-edge 
points are coded as zeros. A listing of a Fortran program that implements 
the edge detection algorithm is included in Appendix B. 

Adjacency Method 

An algorithm called the adjacency method for separation of regions 
of near-homogeneous radiometric data was developed. The algorithm computes 
statistics including the number of observations, sum of data values, 
and sum of squares of the data values for all separate data regions. 

In general, the algorithm uses the following assumptions: 

1. Gradient magnitudes within near-homogeneous data regions (agri- 
cultural fields for these data) are smaller thag near edges 
between regions. 

2. Edges are approximated by gradient magnitudes greater or equal 
to a gradient threshold. 

3. All data points which are smaller than a gradient threshold 
belong to near-homogeneous data regions. 
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4. All data points that are surrounded by an edge and are adjacent 
belong to one and only one data region. 

A flow chart of the procedure is in Fig. 4. In general, the gradient 
magnitude is computed for each data point and a threshold is applied. 

Those data points whose gradient magnitudes are larger than the threshold 
value are regions where edges occur. The adjacency method assigns the 
same numbers to all separate data regions which have gradient magnitudes 
less than the gradient threshold and are surrounded by edges. Statistics 
are computed for each data region using the original radiometric values. 

The gradient threshold is chosen in a similar manner to that for 
the edge detection algorithm. The magnitude of the gradient threshold 
should be sufficiently large to eliminate noise which creates false 
edges but yet retains major edges. The Fortran program which implements 
the adjacency method is included in Appendix B. The resulting statistics 
computed using the adjacency method are analyzed with a classifier de- 
scribed in the next section. 

Statistical Classifier for Field Comparison 

Statistical means and variances for each separable field were 
computed using the adjacency method. Separable fields were those iden- 
tified either by the adjacency method or by an artificial line placed 
in the data where field boundaries were known but were not located using 
the adjacency method. A two-tailed t-test was used to determine which 
field means were statistically separable. A flow chart for the procedure 
is in Fig. 5 and a listing of the Fortran program is in Appendix B. 

Means within a spectral band were ranked and statistically compared. 

Classes of data were generated by combining statistically similar means 
(at 0.05 confidence level) into grand means. Resultant data classes 
were statistically separable using the two-tailed t-test at the 0.05 
significance level. The two-tailed t-test (in contrast to the one-tailed 
t-test) was used to reduce bias because of the initial ranking of the 
means. Orthogonality was maintained because each mean was used only once. 
If the mean was statistically separable from the adjacent mean, it was 
put into a separate class. If there was no statistical difference at 
the 0o05 level, original data values for the mean were combined with data 
values for the adjacent mean resulting in the computation of new statistics 
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Fig. 4 - Adjacency method flow chart. 


















RANK MEANS FROM 
HIGHEST TO LOWEST 



Fig. 5 - Flow diagram of the t-test procedure. 
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for subsequent comparison. 

Example Output 

Figure 6 is an example of the original analog film data from the 
RS-7 scanner of the RB-57F, a gray-level SADE map prepared by quantization 
of the 256 output codes of the digitized film data, an edge detection 
map, and an output from the adjacency method. Field numbers correspond 
to those in Fig. 2a and Table 1. Note that the fields which had no 
apparent emittance variations were not separable with boundary detection 
or adjacency classifier (eg., #12 and 13). An artificial boundary 
was placed between those two fields and statistics were hand computed 
using the values of data listings from the tape. 

SKYLAB S190-A 

The black and white positive transparencies and appropriate 
gray-scale wedge calibration strips from the S190-A were digitized 
into 256 gray levels via SADE. The film/filter combinations for each 
station are in Table 2. 

TABLE 2. FILM FILTER/COMBINATIONS FOR CAMERA STATIONS FOR THE S190-A. 


Position 

Film Type 

S190-A Filter 

Wavelength in ym 

1 

2424 

CC 

0.7^. 8 

2 

2424 

DD 

0.8^0. 9 

3 

So-022 

BB 

0. 6-^0. 7 

4 

So-022 

AA 

0.5^0. 6 


The output of the digitizer operated in log mode was codes ranging 
from 0->-255. These output codes were translated into corresponding film 
densities of the duplicate sensitometry strips. Using procedures described 
by Lockwood (1974), densities were calibrated to the originals with sub- 
sequent calibrations to input radiance levels in cal cm"^ min"^,. The edge 
detection, adjacency method, and t-test algorithms were applied to the 
data to reduce and statistically classify the 30 fields. 



a. RS-7 thermal scanner film, 

Approximate scale 1:32,500 ^ight is cool) 

(light is cool) 



c. Edge detection. d. Adjacency method. 


Fig. 6 - Example outputs of processing algorithms. Field numbers similar 
to those in Fig. 2a and Table 1. 
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SKYLAB S190-B 

The color-infrared film from the S190-B mapping camera was used 
as a base map for analyses of digital products o 

SKYLAB SI 91 

The SI 91 data were acquired over Delta Lake Site 390. Forward look 
angles were not affected by clouds but backward look angles had cloud 
interference. An evaluation of the SI 91 data was conducted to determine 
the angular effects of sensor-measured radiance in the various SI 92 wave- 
bands. The radiance variations included varying thicknesses of the atmo- 
spheric and geometric variations in sensor->-target->-sun angles. 

A plot of original data calibrated to radiance (Juday 1974) versus 
wavelength for the normal -to-horizon (0°) look angle was prepared. In 
addition, variance of the radiance measured at the detector which was 
associated with differing atmospheric path lengths and target- to-Sl 91 
look angles was evaluated. The data take required 1 minute and 9 seconds. 
Sensor drift of =^0,15% per minute did not significantly affect the com- 
parison of radiance measurements within the small observation time; 
however, absolute radiance may have been in error because the actual 
timing of the prepass autocal was not known. The relative atmospheric 
path was determined by geometrical configuration of the SKYLAB workshop 
to the location of the ground targeto The shortest atmospheric path 
was at nadir and was given a zero length. The longest path for the SKYLAB 
altitude was 149 km greater than the assumed zero path. The 149- km length 
was calculated by assuming that atmosphere extended to the altitude of 
the SKYLAB which, of course, is not true. However, the relative radiance 
changes from normal in percent of total length are not affected by this 
assumption. A regression analysis of the relative path length to radiance 
was performed. 

SKYLAB SI 92 

Visual analyses of the DAS products from the SI 92 scanner were 
performed to identify which available spectral regions could be used 
to identify water-related soil variations. Digital analyses were completed 
including application of the edge detection and adjacency method algorithms 
for field location and computation of statistics, calibration of the 
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data to radiance, and correlation and regression of radiance data to ground 
variables. The t-test procedure was used to determine which fields were 
separable within each spectral region. The Channel 15-16 thermal data were 
used in preparation of an evapotranspiration map. Specific assumptions used 
in developing the map will be outlined when the map is presented later 
in this document. The spectral regions available for analyses are presented 
in Table 3. The X/5 detector was operative in the scanner during this 
period and this deleted some of the channels available for data acquisition. 
In addition, the channel 11 data were not suitable due to noise. 

TABLE 3. SPECTRAL REGIONS AVAILABLE FROM THE. SI 92 SCANNER. 


Spectral Region 

Channel Number i n ym 


High or 
Low Rate 


3-4 

0.56-0.61 

high 

7-8 

0.68-0.76 

high 

9-10 

0.78-0.88 

high 

19 

0.98-1 .08 

low 

20 

1 .09-1.19 

low 

12* 

1 .55-1 .75 

high 

13-14 

2.1-2.35 

high 

15-16 

10.2-12.5 

high 


* Channel 11 not suitable for analysis 


SKYLAB SI 93 

No SI 93 data were received for the test site. 

SKYLAB SI 94 

Due to the resolution of S194 and the small test site, the S194 
data were not analyzed. 

RB-57F Aircraft 

Aircraft data acquired during this investigation are listed in Table A-2, 
Appendix A. Photographic products were used for positive ground identi- 
fication and illustration. Thermal data from the RS-7 scanner were visually 
interpreted to aid in identifying within-field variations which were 
below the resolution of the S192 sensor aboard SKYLAB. In addition. 
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data processing algorithms were initially tested on the SADE digital 
product of the RS-7 thermal analog film. Example products were illustrated 
in Fig. 6. 

NC-130B Aircraft 

Screening film from the multispectral scanner was visually and 
digitally evaluated to reinforce interpretations of water-related terrain 
variations observed using the SKYLAB data. The NC-130B and SKYLAB SI 92 
data provided multistage, multispectral data for concommitant analyses. 
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DATA PRESENTATION AND VISUAL INTERPRETATION 

Ground Data 

Soil moisture measurements and land-use descriptions for each 
of the 30 fields were presented in Table 1. A correlation analysis 
between means of ground measurements and of spectral radiance for 
the thirty fields was conducted. The simple linear correlation coef- 
ficients are listed in Table 4. Soil moisture contents at the various 
depths were all highly correlated. The only significant correlations 
of land use and soil moisture were between both percent green vegetation 
and percent fallow with (h-Z cm soil moisture. However, only 21% of 
the variation in 0-^2 cm moisture could be accounted for with percent 
green vegetation and 14% with percent fallow. Several fallow fields 
were being irrigated in preparation for planting small grain. Many 
fields of vegetable crops were mature and were being harvested. These 
two conditions probably decreased correlations. However, the ground 
conditions did provide a test site favorable for determining the utility 
of reflectance and emittance data for assessing actual soil moisture 
variations without the interference of a high correlation between soil 
moisture and vegetation conditions. 

SKYLAB S190-A 

The four separate black and white S190-A images which were digitized 
for statistical analyses are presented in Fig. 7. 

SKYLAB S190-B 

The S190-B color-infrared photography (Fig. 1) was used as a base 
map. Accurate location and mapping of field boundaries were derived 
from these data. 

SKYLAB SI 91 

Radiance versus wavelength was plotted in Fig. 8 for the 0° (nadir) 
look angle. Note that the scale of the vertical axis differs by a 
factor of 100 when comparing the reflective (Fig. 8a) and the thermal 
(Fig. 8b) spectral regions. The approximate spectral regions of S192 
data are noted. The target scene was range vegetation and was not 
ground inspected due to its distance from the intensive test site. 

The SI 91 land region was not within the coverage of the SI 92 sensor. 



TABLE 4. CORRELATION OF GROUND VARIABLES. 




Gravimetric 

Soil Moisture 



Land Use 



0^2 cm 

2->l 0cm 

1 0-^3 0cm 

Composite 

% Green 
Vegetation 


% fallow 

Soil moisture 








0 2cm 

1 .000 







2 10cm 

0.822** 

1.000 






10 30cm 

0.666** 

0.881** 

1.000 





Composite 

0.783** 

0.952** 

0.980** 

1.000 




Land use 








% green 
vegetati on 

0.461* 

0.157 

0.198 

0.228 

1 . 000 



% dry 
debris 

-0.031 

-0.144 

-0.328 

-0.259 

-0.173 

1.000 


% fallow 

-0.373* 

-0.035 

0.055 

-0.017 

-0.737** 

-0.537** 

1.000 

where * is significant at the 

0.05 level and 

** is significant at the 

0.01 level. 

28df 



ro 

CO 



ORIGINAL PAGE IS 
OF POOR ftUALrm 



d . 0. 5 -»■ 0. 6 uHi 

Station A-6 


b. 0.6 -► 0.7 jjm 
Station A-5 


c. 0.7 -► 0.8 um 
Station A-1 


d. 0.8 0.9 ym 

Station A-2 


Negative prints of S190-A data which were digitized to statistically analyze 
the 30 field test sites. Area is similar to that of S190-B presented in Fig. 1 
Arrcft^s point to intensive test sites identified in Figs. 2a and 2b. 


ro 
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RRDIRNCE(X10-M yw cm-2 sr“i 

0.C.'0 0.02 0.03 0.05 0.05 0.08 0.10 0.11 0.13 0.14 



chlorophyll 


i I 

Spectral Region 
in ym 

a. 0.56 + 0.61 

b. 0.68 ->• 0.76 

c. 0.78 0.88 

d. 0.98 1.08 

e. 1.09 ->1.19 

f. 1.55 -V 1.75 
Q. 2.1 ->2.35 


WAVELENGTHS IN MICRONS 


|a| Ibi j»d>-l |->e-^l 


.06 2.27 2.48 


Fig. 8a - reflective 

Plot of radiance versus wavelength as measured by the S191 at the 0° look angle 
The major atmospheric absorption bands and S192 spectral regions are noted. 
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The various water, carbon dioxide, and ozone absorption bands 
are noted in Fig. 8 and are evident in the measured radiance values. 

The chlorophyll 'a' and 'b‘ secondary absorption bands are evident 
and appear appropriately at 0.660 ym and 0.642 ym, respectively. The 
ozone absorption band (9.4-?-9.8 ym) in the thermal region is present. 

SKYLAB SI 92 and Multispectral Data from the NC-130B and RB-57F 

The DAS display products of SKYLAB SI 92 data are presented in 
Fig. 9. The area is similar to that in Fig. 2a with the same field 
number annotations. Screening film from the multispectral scanner 
of the NC-130B is in Fig. 10. The scanner aboard the NC-130B did not 
cover the center pivot irrigation systems; however, the RB-57F mission 
in Fig. 11 did image the center pivots. 

The twelve separate agricultural fields are labeled in Figs. 9->-ll. 

In addition, two center pivot irrigation systems were located and labeled 
"B" in Figs. 9 and 11. Ground data on land use and soil moisture were 
not acquired for the two areas. However, upon detailed inspection of 
the low-altitude NC-130B photographic data, the center pivots were assumed 
to be operative because the irrigation systems were in place. Note 
the color-infrared print in Fig. 2a. The system on the north and east 
apparently includes a strip of actively growing vegetation. Region "A" 
is a seep caused by an excess of transported water either by overland 
flow or by lateral movement of groundwater. 

Answers to specific questions of interest to the water resources 
or agricultural specialist were sought from the multistage, multispectral 
data. A visual interpretation of the films was conducted and results 
are in Table 5. Only the S192 high rate channels, which had DAS products 
available, and comparable aircraft spectral regions were evaluated. 

A specific comparison includes fallow fields #7 and #13. Field 
#7 had been irrigated two days prior to the overflight and was wet 
at the very surface with no noticeable drying. Field #13 was dry. 

Notice that in the reflective spectral regions in Fig. 10 patterns 
occur within Field #7. These were associated both with normal soil 
variations and with variations caused by land leveling. Reflectance 
anomalies associated with soil moisture are not apparent between these 
two fields at wavelengths less than the 1.533->1.62 ym spectral channel. 





Fig. 9 - Photographic representations of 
the digital data acquired using 
the SI 92, multi spectral scanner 
aboard SKYLAB. The conversion 
of the digital tapes to film was 
accomplished using the DAS system 
at JSC, Texas. Field numbers 
correspond to those in Fig. 2a 
and Table 1. A similar area as 
in Fig. 2a is included. Note the 
data from the conical scanner 
have not been rectified; there- 
fore, the round pattern associated 
with the center pivots appears 
elliptical (original in color). 

Ca> 

U> 

e. Channel 16, 10. 2-^1 2. 5 um 









n. 1 .533^1 .62 ym 


Photographic representations of the 24-channel multi spectral scanner data acquired with the 
NASA, NC-130B aircraft at approximately 7620-m AGL (25,000 ft) on mission 258. The field 
numbers correspond to the field numbering in Fig. 2a and in Table 1. Dark is high reflectance 
in a -► 0 and is warm for p u. Field numbers are labeled in "i" and *'t". 
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dark is warm, 10.2->12.5 pm 


Fig. 11 




- Thermal data from the NASA RB-57F aircraft mission 260 
collected at 18,300-m AGL (60,000 ft) with the RS-7 
scanner. Field numbers are the same as in Fig. 2a and 
Table 1. Approximate scale 1:50,000. 


TABLE 5. RECOGNITION OF SELECTED LANDSCAPE FEATURES USING SKYLAB AND CORRESPONDING 
NC-130B DATA. 


Landscape Feature 


SKYLAB Data - wavelength in urn 


Comparable NC- 

■130B Data - 

wavelength 

In urn 


0.68-0.76 

0T75-O.88 


2.10-2.35 

10.2-12.5 

0.70-0.76 

0.82-0.88 

i.533-1.82 

2.3-2.43 

11.0-12.0 

Dry fallow vs. wet fallow 
(fields 6, 13 vs. 7,10) 

no 

no 

yes 

yes 

yes 

no 

no 

yes 

yes 

yes 

“Well- or medium" watered 
vegetation vs. “wet" fallow 
(fields 4. 5 vs. 7. 10) 

no 

yes 

no 

no 

no 

yes 

yes 

no 

no 

no 

"Well -watered" vs. “Dry" 
alfalfa (field 4 vs. 3) 

no 

yes 

yes 

yes 

yes 

yes 

yes 

yes 

yes 

yes 

"Medium-^watered" vs. 
"well -watered" alfalfa 
(field 5 vs. 4) 

no 

yes 

no 

- no 

yes 

no 

no 

no 

no 

yes 

Recognition of extent of 
seep (area A) 

no 

yes 

yes 

yes 

yes 

yes 

yes 

yes 

yes 

yes 

Identification of center pivots 
(area B) 

no 

no 

no 

no 

yes 

NA 

NA 

NA 

• NA 

yesl/ 

Survey of irrigated land 

no 

no 

no 

no 

yes 

no 

no 

no 

no 

yes 

Unique location of surface water 

no 

no 

yes 

yes 

no 

' yes 

yes 

yes 

yes 

no 


— The observations are based on visual analyses of the original color DAS products and should be verified using statistical -analyses of 
the digital data. 


y Center pivots not imaged. 


y Interpretation from the 10,2-f-12.5pm thermal scanner (RS-7) on RB-57F aircraft. 
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In the 1 .533->^1.62 ym and the 2.3 2.43 ym channels, a continuous pattern 

transecting the two fields appears to be related to soil variation 
other than soil moisture. The two fields are definitely contrasting 
in all the thermal bands. Field #10 is a fallow field which was under 
furrow irrigation at the time of the data collection. Again, wavelengths 
greater than 1.533-vl.62 ym were effective for separating the wet versus 
dry fallow with no spectral differences apparent in the shorter reflective 
spectral regions. Similar spectral variations in the SKYLAB data are 
evident. 

SKYLAB $193 
No data. 

SKYLAB SI 94 

No analysis of data. 

PRESENTATION OF STATISTICAL RESULTS 

Classification of Fields Using Digitized $190~A Film 

The correlations of the 30 separate field means for each of the 
S190-A films which were digitally analyzed are in Table 6. Each field 
mean was comprised of no less than 20 individual pixel observations. 

The mean values are presented in Table A-3 in Appendix A. 


TABLE 6. CORRELATIONS OF SPECTRAL REGIONS OF THE SI 90-A FILMS. 


Spectral Region in ym 


0. 5-^0. 6 

0.6^0. 7 

0.7->0.8 

0.8^0. 9 

0.5->-0.6 

1 .000 




0.6^0. 7 

0.869** 

1.000 



0. 7-^0. 8 

-0,282 

-0.346 

1 .000 


0.8->0.9 

-0.448* 

-0.567** 

0.366* 

1.000 


where * is significant at the 0.05 level and ** is significant 
at the 0.01 level. 28df 


The greatest redundancy among bands was between the 0. 5-^-0. 6 ym and 
the 0.6^“0.7 ym bands which had an R^ value of 75%, Other correlations 
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were significant but the coefficient of determination was small. A 
statistical classification of fields is presented in Table 7. For 
example, in Table 7 for the 0.7-»<).8 ym region, ten classes of data 
were statistically separable at the 0.05 significance level for the 
31 sites. Fields #4, 25, and 23 were separable from fields #5 and 
11 from fields #22, 19, 24, 28, etc. No statistical separations among 
fields #4, 25, and 23 were present. 


TABLE 7. RESULTS OF t-TEST USING THE MEAN SADE VALUE FOR FIELDS IN 

THE SKYLAB PROJECT. CLASSES ARE STATISTICALLY DIFFERENT FOR 
EACH SPECTRAL BAND USING A 0.05 CONFIDENCE LEVEL. THE FIELD 
NUMBERS CORRESPOND TO THOSE IN FIG. 2 AND TABLE 1. 


Class 

number 

0.7->-0.8 ym 
Station A-1 
Field(s) 

0.8^0. 9 ym 
Station A-2 
Field(s) 

0.6^0. 7 ym 
Station A-5 
Field(s) 

0.8>0.6 ym 
Station A-6 
Field(s) 

1 

4,25,23 

4,5,21 

23 

8 

2 

5,11 

11 

8 

23,1 

3 

22,19,24,28 

17,30,18,15, 

22,8,3,1,2 

1 ,15 

22,13,12,16 

4 

2,3,1,18,8, 

21 

13,24,31 ,19, 
20,23,10,12, 
7,29,25 

16,28,22,14 
17,21 ,13,12 

7,10,15,14, 

2,29 

5 

29,30 

28,26,6 

31,10,25,27, 

7,29,26,2 

27,21,6,19, 

26,28 

6 

20,13,12,10, 

27,16,14 

16,27 

6 

17,30,25,3, 

18,31,24,20 

7 

7,15,6,26 

14 

30,19,20,24 

4,11 

8 

31 

9 

3 

5 

9 

17 

- 

18 

9 

10 

9 

- 

4,5,11 ,9 

- 


The relationship of spectral data to ground information is presented 
using correlation coefficients in Table 8. The 0.6-f0.7 ym spectral 
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band accounted for the greatest amount of variation (36%) in the 0->-2 cm 
soil moisture measurements, A significant correlation with the 2->-10 cm 
soil moisture was also obtained in the 0.6-)-0.7 ym region; however, 
note from Table 4 that the two depths of soil moisture were highly cor- 
related. Therefore, a correlation with one variable implies a correlation 
with the second. The same is the case for correlation with the composite 
moisture because the composite moisture is composed of the individual 
soil layers combined. The 0.7->0.8 ym region provided highly significant 
correlations with variations of percent green vegetation and percent 
fallow soilo The 0.6^0. 7 ym region of the S190-A sensor would be used 
for soil moisture predictions with the 0<.7-^0.8 ym band for estimates of 
land cover if only these data were available. 

TABLE 8. CORRELATION OF S190-A TO GROUND VARIABLES. 




Spectral Region in ym 



0.5^0. 6 

0.6^0. 7 

0. 7-^0. 8 

0.8->0.9 

Soil moisture 





0 2cm 

0.404* 

0.598** 

-0.322 

-0.411* 

2 10cm 

0.222 

0.371* 

-0.166 

-0.127 

10 30cm 

0.258 

0.296 

-0.142 

-0.047 

Composite 

0.279 

0.368* 

-0.177 

-0.115 

Land Use 





% green 
vegetation 

0.565** 

0.561** 

-0.675** 

-0.367* 

% dry debris 

-0.233 

-0.058 

-0.073 

0.039 

% fallow 

-0.324 

-0.440* 

0.628** 

0.287 

where * is significant at 

the 0.05 level 

and ** is 

significant at 

the 0.01 level . 

28df 





Effect of Look Angle and Atmospheric Path Length on S191-measured Radiance 
Calibrated radiance values from a range of look angles from 45° 
forward to 0° were computed. Radiance values within each equivalent 
SI 92 wavelength interval were averaged. Radiance versus atmospheric 
path length was evaluated. The atmospheric path length at 0° look 
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angle was assumed to be of zero distance because the actual atmospheric 
layering of the various absorbing and scattering constituents was unknown. 
At a 45° look angle for the SKYLAB altitude, the relative increase in 
path length from the sensor to the land surface was 149 km when compared to 
the distance at nadir. Therefore, the regression of changing scene 
radiance, which was computed using the SI 91 sensor as the dependent 
variable and changing path length as the independent variable from 
(H149 km, should be linear with respect to atmospheric thickness of 
the various constituents. Sensor-measured radiance changes with differing 
look angles are not only a function of atmospheric path length but 
are also a function of sensor- to-target-to-sun angles for the reflecting 
scene. Therefore, the regressions include both variance associated 
with atmospheric thickness and geometry. Assuming the presence of 
a Rayleigh atmosphere which had minimal scattering in the longer reflective 
infrared wavelengths and that the geometry variations were not wavelength 
dependent allowed partitioning of the atmospheric and geometric variance 
by subtraction of the slope of the regression for the longer wavelengths 
(the value of 4.5 x 10"^ for 2.12-^2.36 ym range) from the slope of 
the shorter wavebands (see Table 9). If these assumptions hold, the 
remaining value is attributed to the atmospheric component. 

The regressions were significantly different than b^=0 at the 0.01 
level for all wavebands evaluated. The positive slopes for the reflective 
regions or the increasing S191-measured radiance with increasing path 
length indicate the influence of the atmospheric scattering. The negative 
slope of the thermal band reveals the cooling effect of the atmosphere 
when the effective atmosphere is cooler than the radiating surfaceo 
Relative atmospheric radiation was 44 times greater in the 0.41->0.46 ym 
spectral region than in the 1.56 ^-1.74 ym spectral region. Regression 
coefficients progressively were larger from the longer reflective-infrared 
wavelengths to the shorter wavelengths in the visible spectrum except in 
the water absorption region of 1.10->-1.18 ym where absorption probably re- 
duced the measured scattering effect. 

Classification of Fields Using SI 92 Data 

The correlation among spectral regions of the S192 scanner using 
paired observations of the means from each of the 30 separate fields 
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TABLE 9. REGRESSION COEFFICIENTS OF THE SI 91 SENSOR TO ATMOSPHERIC 

THICKNESS AND SENSOR GEOMETRY FOR THE SI 92 SPECTRAL REGIONS. 


Spectral band 
in ym 

Slope of Regression 
As Function of 
both Atmosphere 
and Geometry 
(b) 1/ 

Slope of Regression 
As Function of 
Atmosphere 
(b-0.045) U 

Rel ati ve 
Atmospheric 
Radiation 
/b-0.045\ 

V 0.084 / 



0.41 0.46 

3.721** 

3.676 

43.8 

0.46 -)■ 0.51 

3.283** 

3.238 

38.5 

0.52 0.55 

2.903** 

2.858 

34.0 

0.56 ^ 0.61 

2.340** 

2.295 

27.3 

0,62 0.67 

2.235** 

2.190 

26.1 

0.68 ^ 0.76 

2.200** 

2.155 

25.6 

0.78 ^ 0.88 

1.995** 

1.950 

23.2 

0.98 1 .08 

0.967** 

0.922 

10.9 

O 

CO 

0.428** 

0.383 

4.5 

1 .20 1 .74 

0.498** 

0.453 

5.4 

1 .56 ^ 1.74 

0.129** 

0.084 

1.0 

2.12 2.36 

0.045** 

0.000 


10.2 ^ 12.2 

-0.029** 




where * is for b significantly different than b =0 at the 0.05 
level and ** at the 0.01 level, b is the slope of the regression 
line, Y is the dependent variable (SI 91 -measured radiance), and 
X is the independent variable (relative atmospheric path length). 

ij all exponents are 10"^ 
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are in Table 10. The mean radiance value by spectral region for each 
field is in Table A-4 in Appendix A. Few correlations among various 
spectral regions were significant. The thermal channel of 10.2->12.5 ym 
and reflective infrared 2. 1-^2. 35 ym channel had a coefficient of deter- 
mination of 75%. Visual interpretation of the scanner products revealed 
that both of these channels varied with water- related terrain features 
and would, therefore, probably be correlated with each other. Most of 
the reflective “infrared channels were not significantly correlated. 

A statistical classification of fields based upon the SI 92 channels 
is presented in Table 11. Field groupings within each of the spectral 
bands were significantly different at the 0.05 level using the two-tailed 
t-test. The relationship of the classification to ground data is pre- 
sented using correlation coefficients in Table 12. 

Percent dry debris was significantly correlated with the 1.55->-1.75 urn 
spectral region. The surface 0->2 cm moisture was correlated with the 
two reflective-infrared channels which had wavelengths longer than 
lo55 ym and with the thermal -infrared channel. The near-surface soil 
moisture data of 2-^10 cm and the composite soil moisture data were sig- 
nificantly correlated with the 1.55->1.75 ym spectral region. The thermal 
region significantly correlated with the 10^30 cm and composite soil mois- 
tures. Therefore, to prepare an enhancement map illustrating soil 
moisture variations using the field groups (with associated means and 
standard errors) presented in Table 11, either the groupings for the 
1.55^-1.75 ym or 10.2->12.5 ym spectral region should be included. 

The thirty individual fields were stratified into crop, fallow, 
and range for further correlation analyses. The range category con- 
tained only three fields and that category was dropped. Fallow and 
cropped categories were separated by percent vegetation with the re- 
quirement of greater than 20% green vegetation for the cropped category. 
Fourteen fields were classified as cropped and 13 as fallow. The cor- 
relation of spectral data from the SI 92 sensor to ground variables 
for the cropped and fallow categories are in Table 13. In general, 
the correlations between ground variables and spectral regions for 
the fallow fields were greater than for the cropped fields. For the 
fallow fields, the thermal -infrared region was consistently correlated 
with soil moisture. 



TABLE 10. CORRELATION OF SPECTRAL REGIONS OF THE SI 92 SKYLAB SCANNER USING THE 
30 SEPARATE FIELDS AS OBSERVATIONS. 


Channel 

(s) 




Spectral 

region 

in jjm 





0.56->- 

0.61 

0.68^ 

0.76 

0.78^ 

0.88 

0.98^ 

1.08 

1.09^ 

1.19 

1.55^ 

1.75 

2.1-^ 

2.35 

10.2-> 

12.5 


3-4 

0.56^,61 

1.000 








7-8 

0.68^.76 

0.227 

1.000 







9-10 

0.78->0.88 

-0.155 

0.712** 

1 .000 






19 

0.98^1.08 

-0.378* 

-0.002 

-0.090 

1.000 





20 

1.09^1.19 

0.020 

-0.176 

-0.329 

0.006 

1.000 




12 

1.55->1.75 

0.139 

-0.273 

-0.495** 

-0.151 

-0.129 

1.000 



13-14 

2.1^2.35 

0.868** 

0.115 

-0.245 

-0.270 

0.085 

0.093 

1.000 


15-16 

10.2-^12.5 

0.630** 

0.010 

-0.240 

-0.090 

0.077 

0.027 

0.869** 

1 .000 

where * is 

significant at the 

0.05 level 

and ** is 

significant at the 0.01 level. 28df 
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TABLE n. RESULTS OF THE t-TEST USING MEAN RADIANCE VALUES FOR FIELDS FROM THE S192 DATA. 

CLASSES ARE STATISTICALLY DIFFERENT FOR EACH SPECTRAL BAND USING A 0.05 CONFIDENCE LEVEL. 


Class 

Number 

0.56-v0.61pm 
Channel 3-4 
Field(s) 

0.68^. 76pm 
Channel 7-8 
Fleld(s) 

0.78-1-0. 88pm 
Channel 9-10 
Field(s) 

0.98-1-1 .08pm 
Channel 19 
Field(s) 

1 . 09-1-1 . 1 9pm 
Channel 20 
Field (s) 

1 .55-1-1 .75pm 
Channel 12 
Field(s) 

2. 1-1-2. 35pm 
Channel 13-14 
Field(s) 

10. 2-^12. 5pm 
Channel 15-16 
Fiel d(s) 


1 

28,30,13,12, 
29, 1,20,16, 
27,19, 7, 8, 
24,26 

25,30,17,29, 

16 

25,21,17,14, 

23,30,29,16, 

19,26 

18 

4,16,19, 8, 
26,22,20, 3 

13,29,24,12, 

16 

13,28 

28,13,31,16, 
12,29,27,30, 
8, 3 

2 

31,25, 2,17 

4,23, 8,26, 
19,28,20,31, 
1, 3,21,14 

31,20,28, 3, 
22, 8, 6,27, 
7, 2,11,10, 
4,15, 5,12, 
9,13,18 

4 

1 

1,17,19,25, 
22, 8,30,26, 
3,15 

12,16,30,20, 

29,26 

15,20,26, 1 

3 

10, 5,15,11, 
3,22, 6 

22, 7,27,11, 
15,12, 5, 2, 
13 

** 

17,26,29,22, 

8,27,16,19, 

14 

11 

6,27, 5,31, 
20, 7,11 

27 

24,19, 5,18, 
6 

4 

14,23 

6,10,18 

1 

1,20, 3,15 
2, 7,11,23 
5,31,30 

18,25,17,15, 
12, 5,30, 7, 

28,10, 2, 4, 
23,21 

' 1,31,24,19 
8 

22,14,25 

5 

4 

24 

24 

6,12,10,13 

23,14,31,21, 
2,27,13,11, 
10, 6 

14 

15, 6 

4,17,23,11, 
21 ,2 

6 

18,21 

9 

- 

25,24,21 ,28 

28,24 

18 

3,25,22,17, 

5 

7,10 

7 

9 

- 

- 

9 

9 

9 

11 

9 

8 

- 

- 

- 

- 

- 

- 

7,18,10,14, 

2 

- 

9 


- 

- 

- 

- 

- 

23, 4 

- 

10 

- 

- 

- 

- 

- 

- 

21 • 

- 

11 


- 


- 

- 

- 

9 

- 


cn 



TABLE 12. CORRELATION OF S192 DATA TO GROUND DATA FOR 30 TEST FIELDS. 






Spectral 

Region in 

ym 




0.56-^ 

0.61 

0.68^ 

0.76 

0.7&> 

0.88 

0.98->- 

1.08 

1.09^ 

1.19 

1.55-> 

1.75 

2.U 

2.35 

10.2^ 

12.5 

Soil moisture 
0 ^ 2cm 

-0.116 

-0.251 

-0.360 

-0.034 

0.193 

0.455* 

-0.366* 

-0.479** 

2 ^ 10cm 

0.124 

-0,234 

-0.326 

-0.205 

0.103 

0.533** 

-0.104 

-0,.340 

10 ^ 30cm 

-0.091 

-0.123 

-0.133 

-0.283 

0.064 

0.282 

-0.112 

-0.419* 

Composite 

0.075 

-0.175 

-0.219 

-0.247 

0.093 

0.386* 

-0.145 

-0.426* 

Land Use 
% green 
vegetati on 

-0.359 

-0.033 

-0.031 

0.167 

0.078 

-0.119 

-0.342 

-0.248 

1o dry 
debri s 

0.203 

0.226 

-0.138 

0.009 

0.079 

0.386* 

0.215 

0.277 

% fallow 

0.168 

-0.127 

0.121 

-0.149 

-0.121 

-0.162 

0.145 

0.023 


where * is significant at the 0.05 level and ** is significant at the 0.01 level. 28df 



TABLE 13. CORRELATION OF SI 92 DATA TO GROUND VARIABLES FOR THE 14 CROPPED AND 13 FALLOW FIELDS. 


Spectral Region in ym 



0.56^0.61 

0.68->0.76 

0.7&>0.88 

0.9&>1.08 

1.09^1.19 

1.55-^1.75 

2.U2.35 

10.2-V12.5 

Cropped (12df) 









Soil Moisture 









0 -j- 2cm 

-0.259 

-0.156 

-0.488 

0.030 

0.296 

0.495 

-0.250 

-0.347 

2 ->■ 10cm 

0,186 

-0,026 

-0.369 

-0.307 

0.316 

0.522 

0.102 

-0.131 

10 30cm 

0.303 

0,212 

-0,061 

-0.420 

0.244 

-0.036 

0.120 

-0.233 

Compos i te 

0,216 

0.113 

-0.210 

-0.361 

0.292 

0.18§ 

0.072 

-0.243 

Land Use 









% Green 









Vegetati on 

-0.426 

0.018 

-0.067 

0.288 

-0.117 

-0.208 

-0.375 

-0.505 

Fallow (lldf) 









Soil Moisture 









0 2cm 

0.104 

-0.334 

-0.242 

-0.212 

-0.225 

0.492 

-0.358 

-0.688** 

2 10cm 

0.292 

-0.329 

-0.387 

-0.183 

-0.185 

0.613* 

-0.068 

-0.686** 

10 30cm 

0.127 

-0.319 

-0.271 

-0.237 

-0.207 

0.542 

-0.103 

-0.681* 

Composite 

0.177 

-0.333 

-0.310 

-0.226 

-0.209 

0.574* 

-0.123 

-0.704** 


where * is significant at the 0.05 level and ** is significant at the 0.01 level. 
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Enhancement of SI 92 Digital Data by Statistical Techniques 

Many man-machine interactive procedures allow the interpreter 
to level slice the original digital data to enhance certain scene con- 
trasts which are related to specific thematic interpretations. Equal 
level slicing was presented in the example DAS products in Fig. 9. 

The procedure developed for enhancing the SI 92 data was to base the 
levels for slicing on the statistical classification presented in 
Table 11 and to use unequal levels for slicing. The class means were 
compared with weighting between means based upon the standard deviations 
of the adjacent means. This procedure was used to develop separation 
levels for a variable quantizer program to level slice the digital 
data and to display the data on a color screen via SADE. The ground 
region is similar to that for the DAS products in Fig. 9 and the photo- 
graphy in Fig. 2a, both of which included ground sites #1^-13, The 
color enhanced products are presented for the six different high-rate 
channels in Fig, 12. Since the channel 11 data were not of sufficient 
quality for use, the pixels within a line of data from channel 12 were 
doubled to retain similar geometry as in the other five channels. 

The DAS product for the 10.2-+12.5 ym presented in Fig. 9 had two separate 
levels; whereas, the data in the thermal channel presented in Fig. 12 
has five separate levels (levels determined using statistical t-test 
classification of sites #l-i-13). Based upon the correlation results 
of these channels to the ground variables for the 30 test fields, the 
ground variables most related to the spectral data were: 0.56-^0.61 ym 

for percent green vegetation (r = -0.359), 0.68^.76 ym for 0-^-2 cm 
soil moisture (r = -0.251), 0.78^►0.88 ym for 0->-2 cm soil moisture (r = 
-0.360), 1 .55^1.75 ym for 2->10 cm soil moisture (r = 0.533**), 2.1“+2.35 ym 
for 0->2 cm soil moisture (r = -0,366*), and 10.2-^12.5 ym for 0^2 cm 
soil moisture (r = -0.479**). 

Recognition of Water-Related Terrain Features from 
S190-A and Sl92 Statistical Classification 

The same resource questions were asked of the statistical classi- 
fication of data as were asked of the visual interpretation presented in 
Table 5. The results in Table 14 are based upon observations from 
the same fields as the visual evaluation with the assumption that the 




Fig. 12 - Color-enhanced products of the high-rate S192 channels. 

Levels chosen based upon statistical separation of the 
fields (original in color). 
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TABLE 14. RECOGNITION OF SELECTED LANDSCAPE FEATURES USING SI 92 AND S190-A DATA AND THE 
t-TEST RESULTS OF SCENE CLASSIFICATION AT THE 0.05 SIGNIFICANCE LEVEL INCLUDING 
COMPARISON OF NC-130B AIRCRAFT SCANNER.* 


LANDSCAPE FEATURE SI 92 SPECTRAL R;GION IN vM . SI 90- A SPECTRAL REGION IN vM 


O.SG^.ei 0.6&+0.76 0.7&+0.B8 0.98+1.08 1.09+1.19 1. 55+1.75 2.1+2.35 10.2+12.5 0.5^.6 0.&m.7 0.7+0.8 O.B+0.9 


Dry fallow vs. 
wet fallow 


(6,13 vs 7,10) 

no(no)*** 

np(no) 

no(no) 

■ no(no) 

no 

no(no) 

yes (yes) yes (yes) 

no 

no 

no 

no 

"kell- or mediu»-“ 
watered vegetation 
vs. wet fallow 
(4,5 vs 7,10) 

no (yes) 

no (yes) 

no(yes) 

no(yes) 

no 

no(no) 

no (yes) 

yes (yes) 

yes 

yes 

yes 

yes 

well -watered vs. 
“dry" alfalfa 
(4 vs 3) 

yes (yes) 

no (yes) 

no (yes) 

yes (yes) 

no 

yes (yes) 

yes (yes) yes (yes) 

yes 

yes 

yes 

yes 

medi un-watered vs. 
well -watered alfalfa 
(5 vs 4) 

yes (yes) 

yes(yes) 

no(yes) 

yes (yes) 

yes 

yes(no) 

yes (no) 

yes (yes) 

yes 

no 

yes 

no 

Survey of 

Irrigated** 

Land 

no (no) 

no(no) 

no(no) 

no (no) 

no 

no(no) 

no (no) 

yes (yes) 

no 

no 

no 

no 

Unique location 
of Surface water 

yes 

yes 

no 

yes 

.yes 

yes 

yes 

yes 

yes 

no 

yes 

yes 


* The statistical classes must be separable and when more than one field is Included, the fiel '^ must rank in similar order 
to obtain a yes answer. 

** All preceding answers must be positive before entering a yes. 

*** Closest corresponding NC-130B scanner digital interpretation (0.54+0.58, 0.71+0.76, '.oO+O.SS, 0.98+1 .04, ‘1.53+1. 63, 2.1+2.38, 
and 11.2+11.9 pra respectively. 
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fields were representative of the area analyzed. Alfalfa fields with 
similar canopies and varying soil moisture levels were used for interpre- 
tation of irrigation scheduling. Results similar to those in Table 5 were 
obtained for the discrimination between wet and dry fallow surfaces, 
that is, only the longer reflective- infrared wavelengths and the thermal 
infrared data were effective for spectral recognition. The only data 
which could be effectively used for a survey of irrigated lands were 
in the thermal spectral band. An analysis of corresponding scanner 
data from the NC-130B aircraft is included in Table 14. The statistical 
field classification for the 12 test fields (those in Fig. 2a) is pre- 
sented in Table A-5 of Appendix A. Many additional water-related features 
were separable in the NC-130B data that were not separable in the SKYLAB 
data, especially in the shorter wavelength regions. However, the survey 
of irrigated land could only be conducted with the thermal data for both 
the NC-130B and SKYLAB data. The criteria for a yes answer to this 
question was for all preceding categories to have a yes answer (Table 14). 

Least-Squares Multiple Regression for Prediction of Ground Variables 

An advantage of a multispectral scanner is the automatic registration 
of data from various spectral regions. Therefore, the combined use 
of spectral qualities can be easily attained for prediction of a dependent 
variable from a combination of independent variables of spectral data. 

The field means of spectral information from the spectral regions of 
the SI 92 scanner were used to predict the soil moisture variables. 

Eight spectral regions of the S192 were considered as independent var- 
iables. The computer program used for the analysis chose, as the first 
independent variable, the input which accounted for the greatest reduction 
in error sum of squares* The second variable was automatically chosen 
which, when in combination with the first selected variable, provided 
the greatest additional reduction in error sum of squares. This procedure 
continued until the error sum of squares was either zero valued or 
until all possible independent variables had been selected. The F-test 
in an analysis of variance was then used to determine which independent 
variables should be included in the equation at a given level of signifi- 
cance. The multiple correlation coefficient (R) and coefficient of deter- 
mination (R2) were also program outputs. An example analysis is illustrated 
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in Table 15 for prediction of the composite soil moisture. A summary 
of multiple regression analyses for all 30 fields, the 14 cropped 
fields, and the 13 fallow fields is in Table 16o Only those independent 

variables which were significant at least at the 0.05 confidence interval 

were included in the equation. 

The resulting equation from Table 15 for predicting the composite soil 
soil moisture is: 

y = -10.8 Xi -1.7 X 2 + 0.6 X 3 + 111.9 [17] 

where y = composite soil moisture 

Xi = radiance from 10.2-+-12.5 ym region 

X 2 = radiance from 0.98->-1.08 ym region 

X 3 = radiance from 1.55-vl.75 ym region 

The coefficient of determination (R^) for this equation for the thirty 
observations was 52.4%. 

Equations were developed which accounted for 46-^54% of the variation 
in soil moisture with spectral data when both cropped and fallow fields 
were included in the analysis. The 10,2-^12,5 ym, 1.55-^-!, 75 ym, and 
0.98->-l,08 ym spectral regions were the predominant independent variables 
chosen. For the cropped fields alone, no variables alone (note Table 
12 ) or in combination yielded an equation which was significantly related 
to soil moisture. The spectral regions of 10,2->-12,5 ym, 1,55->1,75 ym, 
and 0,56-^0,61 ym were chosen in that order for prediction of each of 
the soil moisture variables except for the 10->-30 cm depth where only 
the first two were chosen. Coefficients of determination ranged from 
71->-90%. Since no significant equations could be developed for the 
cropped fields and the equations developed for the fallow fields had 
larger correlation coefficients, the fallow fields probably had a con- 
siderable influence on the significance of the equations developed 
for all fields. 

The relationship of the predicted versus actual composite soil 
moisture values using the regression equations is in Fig. 13. The 
actual composite soil moisture was plotted against the predicted soil 
moisture for the fallow fields in Fig. 13a and for all fields in Fig. 13bc 
The actual regressions could not be plotted because three independent 
variables were included in each equation. However, Fig. 13 displays 
the variance in prediction. 



TABLE 15. EXAMPLE ANALYSIS OF VARIANCE OF MULTIPLE REGRESSION FOR PREDICTING THE COMPOSITE SOIL 
MOISTURE WITH SI 92 DATA FOR ALL FIELDS. 

Dependent Variable -- Composite Soil Moisture 


Independent Variable 
Radiance from Given 
SI 92 Spectral Region 
in ym 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean 

Square 

F-Value 

R 

R^(IOO) 

in 

Percent 

10.2-J-12.5 

143.62 

1 

143.62 

11.7** 

* 

18.9 

0.98-V1.08 

164.83 

1 

164.83 

13.5** 


40.7 

1 .5&^1 .75 

89.19 

1 

89.19 

7.28* 

0 . 724 ** 

52.4 

1.09-^1.19 

34.62 

1 

34.62 

2.8 


57.0 

0.68-i^.76 

28.74 

1 

28.74 

2.3 


60.7 

0.56-^0.61 

13.93 

1 

13.93 

1.1 


62.6 

Z.U2.35 

25.93 

1 

25.93 

2.1 


65.9 

0.78^^. 88 

0.79 

1 

0.79 

0.1 


65.9 

Error 

257.29 

21 

12.25 




Total 

758.94 

29 





where * is significant 

at the 0.05 

level and ** 

is significant at 

the 0.01 level. 



cn 
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TABLE 16, SUMMARY OF MULTIPLE REGRESSION ANALYSES FOR PREDICTION OF SOIL MOISTURE CONTENTS 
USING THE SI 92 DATA. ONLY VARIABLES SIGNIFICANT AT LEAST AT THE 0.05 LEVEL ARE 
INCLUDED IN EQUATIONS. 


Land Use 

Number 

of 

Fields 

Dependent 
Variable 
fSoil ^io1sture) 

Independent Variables 

R2/ 

(R^)^ 

bii/ 


b2 

. 

bs 

xa 








No 







. 



. 

Stratification 

30 

0 2cm 

-17.4 

10.2-*12.5 

-2.0 

0.7M.88 



0.684 

45. 8S 



2 10cm 

0.9 

1.55+1.75 1 

-9.5 

10.2+12.5 

-1.3 

0.98+1.08 

0.737 

54. 4S 



10 30cm 

-10.9 

10.2+12.5 ' 

-2.0 

0.93+1.08 



0.659 

43.4X 



Composite 

-10.8 

10.2+12.5 J 

; -1.7 

0.98+1.08 

0.6 

1.55+1. 75 

0.724 

52.4* 

Cropped 

14 

0 -<• 2cm 

none significant { 









2 -*■ 10cm 

none significant i 





* 




10 -*■ 30cm 

none significant j 









Composi te 

none significant 







Fallow 

13 

0 2cm 

-23.0 

10.2+12.5 

1.3 

1.55+1.75 

5.2 

0.56+0.61 

0.921 

84.82 



2 -► 10cm 

-15.9 

lD.2+12.5 1 

1.0 

1.55+1.75 

2.0 

0.56+0.61 

0.947 

89.72 



10 30cm 

-13.5 

10.2+12,5 

1 0.9 

1.55+1.75 



0.843 

/ 1 mZih 



Composite 

-15.6 

10.2+12.5 -1 

i 

0.9 

1.55+1.75 

1.2 

0.55+0.61 

0.906 

82.12 


i/ where bx, are the partial regression coefficients, 

ll where X\, are the independent variables, 

1/ v;here R is the multiple correlation coefficient of resulting equation, and 
h/ where R is the coefficient jf determination of the resulting equation. The equation form is 
y = bixi + bzXz + b^^x^ + c 
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Actual Composite Soil Moisture {% by weight) 



Predicted Composite Soil Moisture (% by weight) 

a. Fallow Fields - 13 Fields 



b. All Land Uses - 30 Fields 


Figure 13 - Predicted versus actual composite soil moisture (0->30 cm) for the test fields 

using the multiple regression eguations with S192 data as independent variables. 
The "o" denotes fallow, denotes range, and denotes cropped surfaces. 
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SOIL THERMAL PROPERTY AND 
EVAPOTRANSPIRATION ASSESSMENTS 

Preparation of Evapotranspi ration Map 

The Jensen-Haise model was used for an estimate of ET, where; 

ET = (0.025 T + 0.08) R [9] 

where ET = potential ET in mm min“^ 

r 

T = air temperature in degrees C 

R = equivalent depth of evaporation of 
incoming solar radiation in mm min"^. 

Incoming solar radiation at the SKYLAB overpass (20:12 GMT) was measured 
with an Eppley pyranometer as 1.020 cal cm"^ min"^ with the air temperature 
as 19,3 C. Therefore, the potential ET prediction is 10.0 x lO-^mm 
min“^. Assume that: 1) the actual ET is equal to the potential ET 

for the field having the coolest surface temperature, 2) the actual 
ET is zero for the field with the warmest temperature, and 3) the 
temperature in the Jensen-Haise model is linearly related to the fourth 
root of the radiance as measured by the S192 scanner. The data set 
was used which covered a similar area as presented in Fig. 2a. Therefore, 
the 13 representative field sites in this region were statistically 
classified using the 10.2^12.5 pm channel. The field groupings were: 

Class I = #13, 2; Class II = #8, 3, 1; Class III = #5, 6: Class IV = 

#4, 11, 2, 7, 10; and Class V = #9. Corresponding class means and 
standard errors of the means were computed with subsequent groupings 
of the SI 92 digital data into the variable quantizer program. The 
Class V data corresponded to field #9 which was a free-water surface. 

The Class IV fields were assumed to have actual ET at the potential 
ET rate as predicted using the Jensen-Haise model (10.0 x 10“3 nim min“^). 
Since the data were acquired under low ET demands and the Class IV 
fields represented high soil moisture regions with water freely available 
for evaporation and transpiration, this assumption was probably valid. 

The Class I fields of #13 and 12 were dry fallow and alfalfa, respec- 
tively, which had considerable soil water deficits. The ET was, there- 
fore, assumed to be 0 mm min“^. The Class I fields corresponded to 
the maximum class temperature with the Class IV fields corresponding 
to minimum class temperature of the agricultural fields. The remaining 
two classes were weighted using the fourth root of thermal radiance 



57 


and had calculated ET values of 8.4 x 10 ~ 3 nui min~i and 5.4 x 10“^ mm 
min“^. The resultant ET map is presented in Fig. 14. This procedure 
was developed to suggest alternative methods for deriving a coefficient 
similar to the crop coefficient which can be used to relate actual 
ET to potential ET. The synoptic, yet detailed, observations from 
satellite altitudes coupled with ET estimates can be useful in irrigation 
scheduling, watershed water budget models, monitoring crop growth and 
yields, etc. A multiple regression prediction using more than one 
spectral region may provide improved estimates since an improved rela- 
tionship to actual soil moisture values for the reported observations 
was attained for equations developed using this technique. 

Prediction of Soil Thermal Properties 

The heat capacities of soils can be derived by summing the heat 
capacities of individual soil constituents per unit volume. DeVries 
(1963) reported that the heat capacity relationship to soil constituents 
could be summarized by: 

L ^ ^ j 

c = 0.46 Xi + 0.60 X 2 + X 3 

where c = heat capacity in cal cm“^ 0 "^, Xi = volume fraction of soil 
minerals, X 2 = volume fraction of soil organic matter, and X 3 = volume 
fraction of soil water. 

A typical volume fraction of soil minerals is 50% and a 4% volume fraction 
of soil organic matter. The volume fractions for a given soil do not 
change appreciably with time; however, the water content does change. 
Therefore, the prediction of changing soil moisture levels can be used to 
predict the changes in soil heat capacity. For fallow fields, a predictive 
soil moisture regression equation using S192 data had a multiple correlation 
coefficient of 0.906. If the soil bulk density and organic matter remain 
constant, the same correlation coefficients apply for the relationship 
of S192 data to soil heat capacity under fallow conditions. Water also 
has a considerable influence on the thermal conductivity of soils (the 
resistance of a soil to the transfer of heat). Therefore, the prediction 
of soil moisture from remotely sensed data provides an estimate of both 
heat capacity and thermal conductivity. 
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Fig. 14 - 


Evapotranspi ration assessment of intensive test site prepared using thei^l radiance 
estimates from the S192 scanner as an input into the Jensen-Haise equation (original 

in color). 


00 



59 


Timing of Thermal Data Collection 

The equivalent blackbody surface temperatures of fallow soils under 
wet and dry conditions are shown in Fig. 15 for one diurnal cycle. The 
wet and dry fields represented are #7 and 13, respectively -- refer to 
Fig. 2a. The precision radiation thermometer (PRT) temperature differen- 
tial between the wet and dry surfaces at 14:12 CST was approximately 6 C, 
well within the thermal resolution of the X/5 S-192 SKYLAB detector. Prior 
to 9:00 CST and after 18:00 CST, the two fields could not be separated by 
surface emittence measurements which emphasizes the importance of the time 
of data collection in surface temperature observations. Note that during 
predawn hours, only a negligible temperature anomaly is present. If the 
energy budget differences were only conductive between the two fields, tem- 
perature differences should have been noted during predawn hours. Daytime 
heating of fallow soil surfaces can serve to indicate a deficit in the 
supply of water to the evaporating surface. Cool evaporating surfaces in 
an otherwise moisture deficient area indicate high soil moisture regions. 

Survey of Irrigated Land With Remote-Sensor Data 

An advantage of remote sensing is the synoptic view which allows data 
collection over large areas in a small time period. For the South Texas 
region in the low rainfall month of January, normally the only high soil 
moisture fields are those which have been recently irrigated. One approach 
to a survey of irrigated lands in this region of ^oy^ rainfall is to identify 
all fields which have been converted from rangeland to cropland. Since most 
of the intensive agriculture in this region must be irrigated, a survey can 
be based upon recognition of the geometric patterns associated with fields. 

A survey such as this can be conducted with data such as the high-altitude 
aircraft photography presented in Fig. 2 or such as the space-altitude photo- 
graphy in Fig. 1. However, many of the fields have very low soil moisture 
because of the seasonal changes in land use and crop growth. A survey of 
irrigated lands using the pattern interpretation of the photographic data 
would, during many seasons, considerably overestimate the land actually 
under irrigation at a point in time. Since the thermal data from both the 
aircraft and SKYLAB altitudes could be effectively used to survey soil mois- 
ture variations, a survey of irrigated lands can be effectively accomplished 
using this spectral region. An example map of the irrigated lands of 



PUT TEMP. (*C) 


Diurnal equivalent blackbody surface temperatures as 
measured with a quantitative hand-held radiation ther 
mometer for wet and dry fallow fields, #7 and #13, 
respectively. 
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the test site is presented in Fig. 16. The light (covol regions) in 
the RS-7 scanner data collected with the RB-57F aircraft include both 
irrigated lands and surface water. In this region, the area covered 
by surface water is small compared to the area of irrigated lands. 

Notice the remaining fields which are used for intensive agriculture 
that are not under irrigation on this date. The area covered by this 
single image is approximately 300 km^ (117 miles^). The same survey 
can be conducted with the SKYLAB data presented in Fig. 17. This single 
image represents a survey of approximately 1000 km^ (390 mile^). 

Use of Remote-Sensor Data for Regional Water Budgets 

Two approaches to monitoring regional soil moisture can be utilized 
with remote-sensor data. One approach is to assess the actual soil 
moisture level by observing crop and soil emittance and reflectance 
variations. The second is to assess evaporative losses as an input 
into the water budget. 

The water budget approach [equation 13] requires an estimate of 
evapotranspiration and soil water flux. Present methods approximate 
the flux with known ground conditions and approximate the uT losses 
using a method such as the Jensen-Haise predictive model. The soil 
water flux term can be considerable for certain soils. Unless the flux 
term can be accurately predicted, effective use of water budgets is 
limited. The ETp models have limitations for predicting actual ET of 
crops under stress due to water deficits but can be used for well -watered 
surfaces where the potential ET effectively predicts actual ET. Thus, 
remotely sensed data can, in conjunction with predictive ET models, 
be used to survey wet and dry regions which may include an irrigation 
district, a watr'shed, or other political or geographic areas. Approx- 
imations of actual evapotranspiration rates on a pixel, field, or regional 
basis may be used to further refine the technique. 

The assumption that the air temperature in the Jensen-Haise model 
could be replaced with an apparent surface emittance, as measured by 
a remote sensor was required. A similar assumption could be incorporated 
for the measurement of R^. The assumptions used to relate emittance 
to evaporative flux using the Jensen-Haise model may be replaced with 
ground control and with measurements of actual evaporative fluxes for 




Fig. 16 - Survey of irrigated lands using the RB-57F, RS-7 scanner data. Area imaged is comparable 
to the photographic image in Fig. 1. The irrigated lands and surface water appear light 
in tone (cool). Approximate scale 1:100,000. 

o» 






Fig. 17 - Survey of irrigated land using a DAS product of the thermal 
band from the SI 92 scanner on SKYLAB. The dark tones are 
' either surface water or’ irrigated lands. Note the arrows 

identifying the location of the intensive test site. Ap- 
ORIGINAL PAGE IS proximate scale 1:200,000 (original in color). 
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equation calibration. The actual quantitative soil moisture prediction 
was suitable on fallow surfaces, but for vegetated surfaces for the 
data analyzed, the accuracy of prediction was not suitable for quanti- 
tative assessment. The use of remotely measured emittance and reflectance 
has promise for assessment of the similar type coefficient as crop 
coefficient which is presently used for relating the actual to predicted 
ET. The repetitive-coverage capabilities of satellite sensors observes 
not only the phenological advances of the crop canopies but also these 
phenological stages at various soil moisture levels. The concept cer- 
tainly merits further investigation. 

Remote-Sensor Data for Irrigation Scheduling 

The use of potential evapotranspi ration prediction models such 
as that by Wiegand and Bartholic (1970) can provide estimates of ET 
for well -watered surfaces. The remote assessments of (net radiation) 
and Tq (surface temperature) are required. If remote sensing were 
used in an irrigation district where the crops were of the same species 
and the same phenological stage, the R^, S (soil heat flux), and T^ 

(air temperature) would be of similar magnitudes. Therefore, the ETp 
would be inversely proportional to a term which is a function of tem- 
perature. With ground control points of known ET , the temperature data 
alone could be calibrated to provide quantitative ET estimates. When 
the land use varies, assumption of equal values for R^ and S would 
not hold. However, as the soil water deficit increases to the level 
where the ET rate is decreased, the use of remotely monitored emittance 
variations for prediction of which fields are under stress may provide 
the input required for irrigation scheduling. 

As soil water depletes, the assumptions required to apply the 
predictive ETp equations do not hold so that the ETp estimate becomes 
invalid as an estimate of actual ET. However, refering to the general 
energy budget equation [6], for a given amount of net radiation on 
surfaces where divergence terms can be neglected, the energy is utilized 
in soil heat flux (S), sensible heat flux to the air (A), and latent 
heat flux (E). As the latent heat flux decreases due to soil water 
deficits, the soil heat flux and sensible heat flux must increase to 
maintain the budget. These tenris are both dependent upon temperature 
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gradients for heat transfer. Therefore, as E decreases with the same 
R^, a larger temperature gradient from the surface must occur to appre- 
ciably increase A and S. This results in increased surface temperatures 
which can be monitored as emittance variations with remote sensors. 
Carson (1961) illustrated that as the {>>45 cm soil moisture decreased 
from 4 cm to 1.9 cm of stored moisture the proportion of R used as 
E (or ) decreased from 8r4>14%. Therefore, greater surface 

temperatures were associated with the dryer conditions for the transport 
of the heat by A and S. The unique feature of heat dissipation by 
latent heat is that the heat is utilized in the heat of vaporization 
with the transfer from the surface as a transfer of vapor. The vapor 
transport is dependent upon vapor pressure deficits in contrast to 
temperature gradients. Therefore, the predictive ET models do not 

r 

require a quantitative prediction of the actual ET rate but can be 
used with ground control in irrigation scheduling with surface emittance 
assessments for determing when soil-water deficits are limiting to 
the crop. 

An example evaluation of utilizing the remote surface emittance 
measurements for irrigation scheduling within a crop species can be 
illustrated using alfalfa fields #3, 4, 5, 11, and 12. Even when sig- 
nificant correlations of SI 92 and soil moisture were not obtained for 
all cropped fields, the alfalfa fields were statistically classed into 
three separate categories — #3, 12 versus #5 versus #4, 11. These 
corresponded with composite moisture values of 15.1%, 18.0%, and 24.1%, 
respectively. Therefore, the emittance variations could be used to 
classify the alfalfa fields into moisture categories. 



Most applications of remote sensing require the relative comparisons 
of ground features for the specific evaluation. For irrigation scheduling 
or soil moisture assessments, the state-of-the-art technology is to 
compare the known with an unknown by using the known to calibrate 
predictive equations (the example multiple regressions in Table 16). 
Therefore, methods of accurately locating and registering ground to 
aerial data and of reducing the aerial data to represent the ground scenes 
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are required. In irrigated regions, the soil moisture variations 
are associated with irrigation scheduling which is based upon field 
variations. Therefore, the method for data reduction developed in 
this activity is well suited to operational programs. The method 
includes location of field boundaries, computing statistics of the 
fields, and using a statistical classification which uses means and 
variances with probability theory for establishing significant differ- 
ences. The use of gradient thresholds to remove high gradient or tran- 
sitional radiometric information from the statistics is especially 
an advantage for comparatively low-resolution satellite data. Data 
should be included from as much of the field as possible to yield an 
optimal representation of the field in question. The pixels at field 
margins represent data mixtures in contrast to observations of only 
one data class. When these points are deleted from the set of data used 
to compute the means and variances, the mean is representative of the 
field without being partially dependent on the adjacent areas. 

An example use of the procedures in irrigation scheduling of similar 
species of crops would include the monitoring of ground control fields 
which are in various irrigation or soil moisture stages. The digital 
data can be quickly reduced to locate fields which are radiometrically 
different than their surroundings, compute the statistics, and statis- 
tically compare the fields with the ground control fields to determine, 
at the confidence level selected by the resource manager, the statistical 
differences of the field means, thereby, classifying the fields as 
to their irrigation requirements. 
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CONCLUSIONS 

1. Wavelengths greater than 2.1 ym were required to spectrally distinguish 
between wet and dry fallow surfaces. 

2. Thermal data provided a better estimate of soil moisture than did 
data from the reflective bands. 

3. Thermal data were dependent on soil moisture but not on the type of 
agricultural land use. 

4. The emittance map, when used in conjunction with existing models, 
did provide an estimate of rvapotranspiration rates. A method to 
use remote thermal data to replace the crop coefficients for relating 
actual ET to predicted ET should be evaluated. 

5. Surveys of areas of high soil moisture can be accomplished with space 
altitude thermal data. If both soil moisture and land use are to be 
surveyed, at least one reflective channel must be included in the 
analysis . 

6. Thermal data will provide a reliable input into irrigation scheduling. 

7. The thermal and spatial resolution of the S-192 X/5 detector SKYLAB 
data is appropriate for monitoring soil moisture and for irrigation 
schedul ing. 

8. Time of data collection for soil moisture surveys should be close 
to midday. 

9. At least the 10.2-+12.5 ym, 1.55->1.75 pm, 0.98->l,08 lam, and 0.56-5-0.61 pm 
spectral regions should be included on operational space-altitude 
scanners for soil moisture applications. 

10. Predictive regression equations using S192 data as independent variables 
had multiple correlation coefficients as large as 0.947 and 0.737 
for predicting soil moisture of fallow and cropped surfaces, respectively. 
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TABLE A-1. SENSOR DESCRIPTION AND PLACEMENT FOR GROUND DATA COLLECTION. 


Sensor 

Measurement 

Depth (in cm)V 

Eppley precision 

pyranometer 

Precision 
radiati on 

0.258-)-2.8 ym 

incoming 

thermometer 

8^14 ym 

surface 

Net radiometer 

total net 

surface 

Mark 16 

radiometer 

0.3M.15 urn 

surface 


0.52-j-l .15 ym 

surface 

Mark IRF 

0.60->-l . 1 5 ym 

surface 

radiometers 

0.71->1.15ym 

surface 

Soil heat flux 

plates 

soil heat flux 

-5 

Soil temperature 

thermocouple 

0, -2.5, -5, -7.5, 
-20, -40, -60, -100 

Air temperature 

shielded thermocouple 

+20, +40, +80, 
+160, +240, +320 

Humidity 

atkins 

+100 


where is below the land surface and "+" is above the land surface. 




TABLE A-2. AERIAL DATA COLLECTION. 







Film/ Filter 







Location 

Platform 

Date 

RC8 

Zeiss 


Hasselblad 



Thermal 

Scanner 

Multi spectral 
Scanner 

RAD/SCAT 






nr L 

2 



5 

6 




Carrizo Springs, 
Texas 

NC-130B 

1/28/74 

So-397 

HF3 

2443 

15 

2402 

58 

2402 

25 

2424 

89B 

So-356 

HF3 

2443 

12 



/ 

/ 

Carrizo Springs, 
Texas 

RB-57F 

1/28/74. 

So-397 

2A 

2443 

510NM 

2443 2402 
12 57 

2402 

25 

2424 

89B 

2424 

87 

2443 

12 


/ 



Timber Lake, 
South Dakota 

HC-130B 

6/13/73 

So-397 

2E 

2443 

510NH 

2402 

58 

2402 

25 

2424 

898 


2443 

12 

So-356 

HF3 

/ 

/ 


Timber Lake, 
South Dakota 

RSI- 

Beachcraft 
Night Mission 

6/9/73 









/ 



Timber Lake, 
South Dakota 

RSI- 

Beachcraft 

6/10/73 



2402 

58 

2402 

25 

2424 

89B 


2443 

15 


/ 



Pierre, 

South Dakota 

NC-130B 

9/18/73 

So-397 

2A 

2443 

12 





2443 

12 

So-356 

/ 

/ 


Pierre, 

South Dakota 

RSI 

Beachcraft 
Night Mission 

9/18/73 









/ 



Pierre, 

South Dakota 

RSI- 

Beachcraft 

9/18/73 



. 2402 

58 

2402 

25 

2424 

898 


2443 

15 


/ 




ro 
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TABLE A-3. RADIOMETRIC VALUES FOR EACH OF THE 31 TEST SITES FROM THE 
SI 90- A CAMERA SYSTEM. FIELD NUMBERS CORRESPOND TO THOSE 
IN FIG. 2 AND TABLE 1 . 


Field No. 

Radiance Values for Each Spectral Region* 

0.5-0.6ym 

0.6-0.7ym 

0.7-0.8um 

0. 8-0. 9ym 

1 

0.125 

0.134 

8.692 

6.462 

2 

0.144 

0.155 

7.934 

6.462 

3 

0.161 

0.172 

8.692 

6.389 

4 

0.176 

0.213 

6.022 

4.472 

5 

0.183 

0.222 

6.551 

4.648 

6 

0.150 

0.155 

11.890 

7.528 

7 

0.140 

0.151 

11.730 

7.063 

8 

0.115 

0.124 

9.098 

6.302 

9 

0.208 

0.216 

24.820 

11.150 

10 

0.142 

0.151 

10.870 

6.967 

11 

0.179 

0.222 

6.641 

4.886 

12 

0.130 

0.146 

10.730 

6.967 

13 

0.135 

0.146 

10.730 

6.626 

14 

0.144 

0.142 

11.540 

8.975 

15 

0.142 

0.137 

11.890 

6.160 

16 

0.130 

0.139 

11.380 

7.916 

17 

0.146 

0.144 

16.350 

5.860 

18 

0.162 

0.176 

8.975 

5.940 

19 

0.152 

0.162 

7.511 

6.794 

20 

0.162 

0.162 

10.410 

6.967 

21 

0.150 

0.144 

9.245 

4.712 

22 

0.134 

0.142 

7.308 

6.231 

23 

0.125 

0.103 

6.317 

6.967 

24 

0.162 

0.164 

7.614 

6.717 

25 

0.159 

0.151 

6.091 

7.143 

26 

0.152 

0.153 

12.080 

7.426 

27 

0.148 

0.151 

11.050 

8.006 

28 

0.152 

0.141 

7.614 

7.426 

29 

0.146 

0.153 

9.655 

7.143 

30 

0.155 

0.162 

9.811 

5.940 

31 

0.162 

0.149 

12.450 

6.794 


* All values reported have exponent of 10"^ for cal cm"^ min"^. 
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TABLE 

A-4. RADIOMETRIC MEAN VALUES FOR EACH OF THE 31 TEST SITES FROM THE S192 
FIELD NUMBERS CORRESPOND TO THOSE IN FIGURE 2 AND TABLE 1. 

SCANNER. 

Field No. 




Radiance /a’ Jes for Each Spectral 

Region* 





0.56-*0.61jiin 
Channel 3-4 

0.6&T0,76um 
Channel 7-8 

0.78+0.88um 
Channel 9-10 

0.98+1 .03 jrn 
Cnannel 19 

1 .09+1 .1 9um 
Channel 20 

1.55+1.75pii: 
Channel 12 

2. 1+2. 35pm 
Channel 13-14 

10. 2+1 2. 5pm 
Channel 15-16 

10. 2+12. 5pm 
Channel 21 

1 

5.181 

5.104 

1.670 

3.558 

2.054 


8. 

295 

0.404 

0.857 


0.856 

2 

4.07S 

3.690 

3.747 

3.217 

1 .347 


5. 

063 

0.150 

0.S01 


0.803 

3 

3.533 

5.004 

4.277 

3.433 

2.214 


7. 

803 

0.285 

0.865 


0.855 

4 

2.587 

5.877 

3.477 

5.841 

2.762 


5. 

053 

0.084 

0.819 


0.818 

5 

3.667 

3.820 

3.229 

2.958 

1.761 • 


6. 

634 

0.261 

0.846 


0.838 

6 

3.413 

3.303 

3.928 

2.539 

1.156 


6. 

950 

0.326 

0.840 


0.837 

7 

4.788 

4.345 

3.796 

3.208 

1.669 


6. 

143 

0.180 

0.796 


0.800 

8 

4.520 

5.625 

4.133 

4.124 

2.338 


8. 

163 ■ 

0.352 

0.867 


0.849 

8 

0.249 ■ 

0.083 

3.183 

0.053 

0.043 


0. 

032 

0.009 

0.749 


0.757 

10 , 

3.790 

2.832 

3.580 

2.366 

1.231 


5. 

304 

0.159 

0.795 


0.810 

n 

3.501 

4.146 

3.641 

3.008 

1.285 


6. 

043 

0.224 

0.808 


0.812 

12 

5«2;5C 

3.844 

3.191 

2.518 

1.777 


8. 

615 

0.470 

0.8S6 


0.874 

13 

5.500 

3.428 

2.691 

1.993 

1.308 


9. 

043 

0.555 

0.892 


0.882 

14 

3.187 

4.974 

8.159 

3.596 

1.590 


0. 

769 

0.156 

0.834 


0.822 

15 

3.654 

4.094 • 

3.398 

3.277 

1.819 


1. 

473 

0.327 

0.859 


0.855 

16 

5.135 

5.918 

5.489 

3.590 

2.494 


1. 

725 

0.465 

0.887 


0.910 

17 

4.061 

6.349 

3.994 

4.523 

1.922 


1. 

671 . 

0.268 

0.812 


0.855 

12 

2.352 

2.244 

2.684 

8.014 

2.047 


0. 

277 

0.163 

0.844 


0.793 

19 

4.323 

5.387 

5.419 

3.601 

2.412 


1. 

655 

0.359 

0.849 


0.737 

20 

5.165 

5.166 

5.097 

3.447 

2.274 


1. 

264 

0.443 

0.858 


0.842 

21 

2.244 

4.999 

9.295 

2.324 

1.466 


0. 

933 

0.017 

0.802 


0.846 

22 

3.533 

4.4S2 

4.170 

4.252 

2.274 


1. 

647 

0.269 

0.835 


0.867 

23 

2.843 

5.770 

7.179 

2.975 

1.644 


0. 

985 

0.116 

0.811 


0.812 

24 

4.517 

2.216 

0.305 

1.310 

6.775 


1. 

783 

0.385 

0.854 


0.858 

25 

4.354 

7.054 

10.132 

■ 1 .355 

1.939 


1. 

652 

0.271 

0.824’ 


0.879 

26 

4.483 

5.488 

5.266 

4.514 

2.296 


1. 

612 

0.431 

0.855 


0.822 

27 

5.059 

4.289 

3.919 

4.062 

1.336 


1. 

345 

0.404 

0.873 


0.826 

22 

5.72? 

5.362 

4.777 

0.987 

0.695 


1. 

106 

0.523 

0.892 


0.867 

29 

5.235 

6.288 

6.322 

4.396 

2.047 


1 . 

814 

0.440 

0.378 


0.841 

30 

5.624 

6.478 

6.567 

2.708 

1.750 


1. 

613 

0.458 

0.872 


0.ES6 

31 

4.380 

5.151 

5.103 

2.804 

1.577 


1. 

280 

0.389 

0.889 


0.863 


*AH Values reported have exponent of 10'^ for cal cm-^ ntin"^. 
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TABLE A-5 RESULTS OF THE t-TEST USING MEAN RADIANCE VALUES FOR FIELDS FROM THE SCANNER OF THE 
NC-130B AIRCRAFT. CLASSES ARE STATISTICALLY DIFFERENT FOR EACH SPECTRAL BAND USING 
A 0.05 CONFIDENCE LEVEL. 


CV.zz 

JiLT-.er 

Ch-1 . 
FieHfs) 

. 0.4U 
O.AS-^rt 
Ch-2 
Fleldfsl 

0.45* 

0.S2;*!I5 

Cfi-3 

Fleldfs) 

0.54-^ 

O.S8ua 

C.h-4 

nel3(i) 

C.59- 

0.64bm 

Ch«5 

Fleld(s) 

0.65^ 

0.69un 

Ch-6 

nsii(5i 

0.7U 

0.76'jtB 

Ch^7 

Ffeldfs) 

0.77- 

0.31*.im 

Ch-8 

Fleldfsl 

0.83- 0.»- 1.20* 

0.88utn 1.34t9 1.30i;m 

Ch-9 Ch-lC Ch*n 

fleldls) Fie!d(.-) Field(s) 

1.53* 

l,53tJflt 

Ch-12 

Ffeldfs) 

2.10- 

2.38Mn 

Ch-13 

Fleldfsl 

8.5- 

3.9uni 

Ch-17 

Fleldfsl 

9.5- 

10.2ufi 

Ch-19 

Ffeldfs) 

10.2- 

ll.Oym 

Ch-20 

Ffeldfs) 

11.2- 

ll.Sim 

C.V21 

Fleldfsl 

12,2- 

13.0:ja 

Ch-22 

Ffeldfs) 

1.14* 
1.1 5vn 
Ch-23 
Ffeldfs) 

1.05* 

T.C?yK 

Ch-24 

Ffeldfs) 

1 

S 

7. 3 

8. 7 

7 

7. 8 

7. 8 

4,11 

1.11 

11, A 

11. 4 

4 

1 

8 

12 

12 

13 

13,12 

13.12 

1.11 

5. 4 

2 

1 

1.13 

1.13 

e 

1.13 

1,13 

5 

5 

5 

e 

0 

6 

13; 1 


3,13 

12 

3. a 

8 

5.1.7.8 

n. 5 

3 

10 

6 

6 

1,13 

6. 2. ■ 
10 

2, 

12 

7. 2. 
3 

3 

3. 7, 

1 

7. ^ 

11. 5, 
7 

13 

6.12 

13 

1. a 

3, 8 

1 

1. 3 

13. 3 

7 

4 

13. 6 

10.13 

10,12 

6. 2 

u 

10 

fi. 1 

7. 2, 
1. 8, 
13.12, 
6,10 

2. e 

e 

13. 1 

12 

3 

1 

6. 5 

1 

6 

6 

12,10, 

2 

s 

5 

Z. 7 

3 

3 

10,12 

3 

3 

13.12 

- 

13,12 

1. 2. 
13 

3 

6. 3 

I. 2 

8 

4, 2 

6 

5 

s 

6 

3.12 

6 

n 

3 

3 

3 

4 

4, 5 

10. 6 

- 

10 

U 

12, 2 

2. «. 
7,10,5 

10 

6. 5, 

7.1,0. 

n 

5 

4 

4 

. - 

13. 2. 
10. 1 ■ 

7 

3. < 

■ 4' 

4 

4 

5 

11 


- 

6 

10, ( 

10, 6 

11 

4. 5 

4 . 

- 

4 

2.11 

11. 2 

- 

- 

Z 

11. S 

s. s 

s 

S 

11 

- 


- 

- 

- 

- 

- 

11 

11 


n .-2 

10 
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- 
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9 

- 
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11 
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EDGE DETECTION COMPUTER PROGRAM 

The subroutines required to implement the edge detection program are 
MASMOV, HISPLT, and SYMFRQ, which may be obtained upon request. 


C PARAMATER CARD 

DII^ENSIQN IRBI2) 

INTEGER*^ LiNEPK(5l2),IOFREQl256),IGFReQ(256) 

INTEGER+A ALABEH20),MEDG4(512) 

INTEGER<=2 M2 
LOGICAL=!‘l LINE,LIN3( 2) 

LOGICAI^l LINK 2048), LIN2( 2048) ,MEDG1( 2048) 
EQUIVALENCE ( MEDGE ( I ) , MEDGl I I ) , MEDG4 1 1 ) ) 

INTEGER+2 MU N I ( 1024) , ML IN2 1 1024 ) t ML IN3 U ) 

EQUIVALENCE ( ML IN3 ( 1 ) , L IN3 < 1 ) ) 

DIMENSION NUMK 16) ,NUM2( 16),NUM3( 16) 

INTEGER<<2 I D I R { 1024, 10 ) , MRBGR2 ( 1024, 10 ) , MEDGE (1024) 
INTEGER«2 NSUM2(1024) 

INTECER>!‘2 NSUMK 1024) 

COMMON /KLENG/ILEN/KCALL/ICALL 
COMMON / INPUT/LINE (2048) ‘ 
equivalence (HUNK 1) ,LIN1( 1) ) 

EQUIVALENCE ( ML IN2 ( 1 ) , L IN2 ( I ) ) 

EQUIVALENCE ( L I NE ( 1 ) , L INEPK { 1 ) ) 

DATA NUMl/1, 2, 3, 4, 5, 1,2, 3, 4, 5, 1,2, 3, 4, 5,1/ 

DATA NUM2/ 2, 3, 4, 5, 1,2, 3, 4, 5, 1,2, 3, 4, 5, 1, 2/ 

DATA N'.!MT/3,4,5,l,2,3,4,5,l,2,3,4,5,l,2,3/ 

DATA rnFRFQ/256*0/, !GFREQ/256«0/ 

ICR=11 
LP-12 
MAG I N= 14 
MAGOU-15 


C .MXI ...STARTING HORIZONTAL BYTE NUMBER. 

C..........M.XF.... .STOPPING HORIZONTAL BYTE NUMBER. 

C ....MYl STARTING VERTICAL LINE NUMBER. 

C. ...... ...MYF STOPPING VERTICAL LINE NUMBER. 

C. .NFRQ. .. .CONTROL FOR HISTOGRAM OUTPUT. 

C.......... 1...AND...2 NO INDIVIDUAL COUNT ONLY BAR GRAPH. 

C.... 3. ..ONLY INDIVIDUAL AND CUMLITIVE COUNT NO BAR GRAPH. 

C.......... 4... AND ABOVE BOTH BAR GRAPH AND COUNT. 

C. ..LRBGR.. .GRADIANT THREHOLO . THAT WANT TO KEEP PERMANNATLY ON 

C.. FILE, ANYTHING LOWER THAN THIS IS SET TO ZEROIO). 

C..... ... ..NPRT,...IF GRADIANT OUTPUT IS NOT WANTED SET EQUAL TO ONE, 


C.......... . OTHERWISE ANY NUMBER INCLUDING ZERO. 

READ ( ICR ',1000) MXI ,MXF, MY I, MYF, NFRQ, LRBGR, NPRT,NDIRGR 

MXF4=MXF/4 

MXF=MXF4*4 

NR=MYF-MYI 

GO T0(l08) ,MY1 > 

MY1=MYI-1 • 

DO 100 I 1=1, MY I 
' 100 READIMAGiN.lOOl) 

108 1LEN=MXF-MXH-1 • 

KM=ILEN-1 • 


OEIOINAB PAGE 18 
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KML=ILEN-2 
NNN=3l 
ICALL=l 
IYY=MY I+,2. 

MX=MXl+l‘ 

NST=NXI-l . 

ILEN2=ILEN><‘2 
K4=KM/4 
111 = 0 

DO 104 1=1, NR 
III = MI + 1 

READ (MAG IN, 1001) ( L I NEPK ( J ) , J = I , MXf4 ) 

CALL MASMOV( ILEN,LINE(MXI ) ,LINE( 1) ) 

LIN3(2)=LINE(1) 

C., SUMS TWO PIXELS TOGETHER ACROSS THE LINE. 

DO 51 1T2=4, ILEN2,2 

IT=IT2-2 

m'=IT/2 

ITR=ITT+1 

M2=MLIN3a) 

LIN3(2)=LINE(ITB) 

51 NSUM2 ( ITT)=M2+ML IN3( 1) 

GO TO ( 14,15, 16, 10, 10,10,12,12,12,12,12,12.13) ,I 1 1 
13 I I 1 = 8 

C SECTION FOR LOCATING EDGES 

12 IBND=III-l 
I1 = NUM1( IBND) 

I2-NUH2{ IBND) 

I3 = NU.‘13( 13ND) ■ ' 

IYY=IYY+l 

C. CHECK EACH POINT FOR LOCAL MAXIMUM ACROSS THE LINE. 

DO 22 IT=2,KKL 
ICX=0 

ITPLUS= IT+1 
ITMIN=IT-1 
IDX=I0IR (IT, 12} 

IR6M=MRBGR2( IT, 12) 

'IF{ IRBM.LT.LRBGR)GQ TO 25 
GOTO {24,28,26,23,24,28,26,23), IDX 
C, ....CHECK FOR LOCAL MAXIMUM IN DIAGQNNAL DIRECTION. UR TO LL 

23 IRB1=IRBM-MRDGR2I ITPLUS,I1) 

IRB2=IRBM-MRBGR2( ITMIN, 13) 

IF(lRBl) 22,41,41 

41 IF( IRB2) 22,61,61 

61 CONTINUE. 

GO TO 21 

C. ....CHECK FOR LOCAL MAXIMUM . ABOVE AND BELOW 

24 IRni=IRBM-HRBGR2( IT, ID 
IR82=IRBM-MRBGR2( IT, 13) 

IF( I RB I ) 22,43,43 ’ 

43 IF( IRB2)22,62,62 ' 

62 CONTINUE ’ 

GO TO 21 

Ci.. ...... .CHECK FOR LOCAL MAXIMUM DIAGQNNAL DIRECTION . UL TO LR 
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28 IRBl=IRBM-MRBGR2nTMIN,ll) 

IRB2=IRBM-MRBGR2( ITPLUSi 13) 

IFURBl) 22,45,45 
45 IF( IRB2)‘22,63,63 

63 CONTINUE 
GO TC 21 

C **CHECK FOR LOCAL MAXIMUM. RIGHT AND LEFT 

26 IRB1=IRDM-MRBGR2(ITPLUS, 12) 
iRB2=IRBM-MR8GR2( ITMIN, 12) 

IF( IRBl) 22,47,47 
47 IF( IRBZ) 22,64,64 

64 CONTINUE 

21 ICX=IRBM 

C MAXIMUM GRADIANT VALUE IS 255. 

IF( ICX.GT.255) ICX=255 
■ 25 CONTINUE 
2 CONTINUE 

22 MEDGEl ITM1N) = ICX 
DO 73 IP-^2,ILEN2,2 

73 MEDGK IP/2) = MEDG1(IP) 

MX=l 

NNN=31 

GO T0(986),NPRT 


C ....WRITE OUT GRADIANT TO TAPE. 

WRITEtMAGQU, 1001 ) {MEDG4I IT) , IT=1,K4) 

986 CONTINUE 

CALL SYM FPQ ! 0, ITMIN s MFDQF n. ICFREQ) 

C GRAPI FNT SFCTION “ ’ . 

C DETERMINE DIRECTION AND GRADIANT FOR EACH POINT. 


10 Il=NUMl( III-3) 

DO 52 IT=l,KM 
IT2=IT+1 

X=M1.IN1( IT2)-MLIN2IIT) 

Y=MLIN1 ( ITl-ML IN2( IT2) 

IDIRI IT, U)=IDIRX(X,Y) 

GGRAD= ABS(X)+A8S(Y) 

■ 52' MRBGR2I IT, Il)=GGRAD 

C..., MOVES AVERAGES FORM ARRAY TWO TO ARRAY ONE. 

16 CALL MASMOVI ILEN2,L1N2I 1) ,L1NU in ' 

C...., ADD SUM OF TWO ROWS TOGETHER DIVIDE BY FOUR. AVERAGES .FOR A 

C........ .. FOUR PIXEL ARRAY, 

15 DO 74 IT=1, ILEN ■ 

74 MLIN2I IT)=INSUM2( IT)+NSUMl(IT))/4 

C.... MOVES SUM FROM ARRAY TWO TO ARRAY ONE. 

14 CALL MASMOVI ILEN2,NSUM2( 1) ,NSUMl( 1) ) 

104 CONT INUE 

WRITe(LP,i009)ITMIN, lYY 
IGFREQI l) = 0 

READ(TCR,1003,END=110) ALABEL 
C,. ....... .PLOTS HISTOGRAM OF OUTPUT. 

110 CALL HISP.LTI IGFREO, ALABEL, NFRQ) 

'ICOO F0RMATI815) 

1001 FORMAT (2 00 1 1 2A4) ) . 

1002 FORMAT! ' • ,3214) 
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1003 FQRMAT(20A4) 

1004 FORMATC ’ 1* I l0Xt20A4) 

I0n«j FORMAT(‘ » ,5X, 'ELEMENT ',15, 'LINE », 1 5 , 73X, • ELEMENT 15, 'LINE 

15) 

1006 FORMAT {200(50A1 ) ) 

1009 FORMAT!' ',10X, 'THERE ARE', 14,' ELEMENTS USED, AND ',14,' LINES') 
END 

FUNCTION IDIRX(X,Y) 
iF{.NQT.(X.EQ.O.AND.Y.EQ.O))60 TO 5 
IDIRX=0 
RETURN 

5 IF(X.EQ.0)X=.001 
RAT=Y/X 

IFI.NOT. (RAT.GT..4142.AND.RAT.LT.2.4142) )GQ TO I 
IDIRX=2 

IF(Y.LT.O) IDIRX=6 
RETURN 

1 IF( .NOT. (RAT.LT.-.4142.AND.RAT.GT.-2.4142) )60 TO 2 
IDIRX=4 

IFIX. GT.O) IDIRX=8 
RETURN 

2 RAT=ABS(RAT) 

1F(RAT.LT.2.4142)G0 TO 3 
I0IRX=3 

IFIY. LT.0) I0IRX=7 
RETURN 

3 !DIRX“1 

IF(X.LT.O) 1DIRX=5 ■ 

RETURN 

END 
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ADJACENCY METHOD COMPUTER PROGRAM 

The subroutines required to Implement the adjacency classifier program 
are REDAPE, MASMOV, CLSAPE, which may be obtained upon request. 

*************** t***^****!/*************************************** 

C I’HCI.KAH NAUR... ADJACtNCY CLASS IT lEH 

i: DcrniuM 2,1 ’>as 

C. WMinUI Hi MIKE KUSSLLL 

************ *4* t*$»**********0****ir¥******4*******n***** 

C CCNThCL StCMON 

ntl-'t.MS lull MADJI5121 
OIMCNSlfJN LINSUMI lOZ'V.aj 
(Nr(;OLH*A L$1256) 

INtl:CLK»'', MUH,SIJt',SUH2 

tNrEGLR«2LI'JErL|Nll 1024) »LA»DEC,0ECt(102<V) tllNOUPI 1024) 

LOGICAL*! LC,A(bl2) 

COMMON NUM(256.2)iSUM(2b6,2),SUH2(256r2l, IBG.IC0iII,LA(2048) 

COMMON OLCi 1024) iLCI 256) 
common /I'lPur/LINEI 1024) 

EOUIVALCNCE { DSC 1< I ) ,LC « I ) ) 

EOUIVALLNCC (Lfiin.Am) 

EQtJlVALKf '*fi ILINCUPt n.LINEin) 

C INiriAtlZE CONSTANTS AND SET LIMITS 

r,********************************************************************* 

C OFFSET INTO TRANSLATION TABLE 

ITAD=0 
C 

C ICO IS THE STARTING CLASS 

C CLASSES 1-4 ARE RESERVED FOR EDGES 

C CLASS 5 IS RESERVE) FOR INSIGNIFICANT CLASSES 

tCD=3 
inG=ICD ■ 

C 10 S IDO ARE VARIABLES TO CONTROL SEARCHING TECHNIOUECIFLC) 

1 0=1 
100 = 2 

C ICR=m)MDER OF CARD READER 

C LF= NUMBER OF WRITER 

C MDIN=NUM3ER OF TEMPORARY FIELD NUMBER STORAGE AREA 

C MDOUT=NUMBER OF PERMANENT FIELD NUMBER STORAGE AREA. 

ICR = S 
LP=6 
M0IN=8 
MOOur=io 

C IERR=l IMPLIES TRANSLATION TABLE IS FULL AFTER TRANSLATION 

IERR=2 

C INITIALIZE TRANSLATION 

LA11)=1 
LA(21=2 
LA(3)=3 
LA14)=4 
LA15)=5 

1004 FORMAn'O',' LINE*'/32I4I 

1005 FQRHATCO',* IFLG= ',110,' NUTS=»,I10) , 

C READ INITIAL ‘DATA BLOCK LIMITS I6I5I 

READ! ICR aCOQIHXl.HXF.MYI .MYF.LRBGR.IAVE 
C H>;i,MXF IS the LUCATIUN UF I HE INITIAL POINT 

c rvr is the lc'cayio*' of the final point • ‘ 

C LROCa IS THE GRADIENT THRESHOLD 

C NREC IS THE NUMBER OF ROWS MINUS ONE FROM THE DATA 

C lAVE IS THE NUMBER OF DATA LINES AVERAGED 

NREC=MYF-MYI 

C ZERO THE NUMfSUH.e SUM2 

DO 200 1=1,256 
NUHl 1,1)=0 
SUM(I,lJ=0 
200 SUN2! I ,11=0 
C . 

G START THE LOOP OF READING DATA 

C******************************************************************^** 

C SKIP PYl DATA RECORDS AND READ ONE RECORD INTO CQMNON/TNPUT/ 

DO ICO S l = l,MYI 
100 CALL REOAPEIMXI.HXF.KHT ' 

K=KM ’ 

MZ=2*K 

IIELEM=IAV6*IAVE 
11=0 



82 


C Jveiuce tAVE by IAVE .MEIGHOORHODlIS FCK first UNf:,, 

GO TG( 52 )tlAVE 
I 1 AV=IAYE -1 
K'^KM-i AVE 
fU=2*K 

00 ‘51 1 = 1 , IAVE 
no BO J=l,K 
JMAX=JMIAV 
no 50 KK=J,JHAX 

50 UNSUF(.l,n = UNSUMl J, IlHlMElKKl 

51 CALL KEOAP£(FXI,FXF,-KHI 
00 55 J= 1 ,K 

SUS = 0 . 

DO 5 A l«l.lAve 

5 A SUS = SUS^LINSU«( J ,1 1 
LIM UT=SUS/I lELEM 
55 COiiriNm 
GO TC 53 

52 CAIL F'ASWVIME.lINOUIM. INI) 

5 J M) =N- I 

KI' = K -1 
N/. = 2*M'. 

1 1 ^0 

C... ...WPfN TONE WITH OATA 

00 HI UiHAVIHNKhC.lAVt 
00 101 KNHvlilAvF 
n-t 1 t 1 

|,<UI 111 I’Al'l I.MXIvMiirtKM) 

r.... ...AVI lAGt Ni M I INI 

r... 4 ..im OnA.NCH AKOHMI I INI. AVIllAGF Sll'SltlN 


Gil n:iH/ii,iAvn 
|)fl 5/ .lj4l,K 
I.IM'-UFU J,KKUl =0 
JHAX = JJ+nAV 
UO 5il KK=JJ,.IMAX 

5(1 LIN'jUfMJJ,KKHI«L INSOMnJrKKlD+Lir.eiKK) 

SlJS=-0 

L'lJ 50 KB’IiIAVE 
05 yj5 = 5U'>«UNSUM(JJ,Kni ■ 

5/ Llfll ( J.J1=SUS/I1CL£H 
773 IFLG=IO 
NAIJ.J = A 
IFL= 1 
ICI1NU('«2 

C.4. ...OFOCCSS A LIHE OF DATA FOR ADJACENT CLASSES 
DO 102 J=l,KO 
J2=J+1 
N=UNUJI 
Y=N-LIM£(J2I 
X = Um IJ2J-LINEI J) 

C APPROXIAHATE KAGNITUDE OF GRADIENT VECTOR.. 

Z=(XA>X»Y'»Y)F’».5 

iKBCR=Z 

C... ...CCMPAHE GRAOfENT MAGNITUDE TO THRESHOLD.... 


IF( IRHCR-LRDGft) llO, 111, HI 

C... ...IF GRADIENT IS LESS THAN THRESHOLD....... 

110 CO TUI U3,U4,U5,U5>,IFLC 

C.,. ...SEARCH r-OSiriONG A and C FOR .....FIELDS C... 

C. FIELD NUMBERS ADJACENT TO B A 8... 

C SECTION IFLG=3 OR 4 , 

U5 IB=0£CU J) 


C... ...CHECK POSITION C FOR INSIGNIFICANT FIELD NUMBERS....... 

GO rciSO'J, 505, 509,509, 516), IB 

C YES THERE IS AN ADJ CLASS (CHECK IF SAME AS BEFORE) 

GO rc( 506) , IFL 

C CONTINUE INC NEW ADJ) 

GO TC 114 

516 GO 10(519,519,519, 519, 519), NAOJ 
GO TO 114 

519 NADJ=5 

GO TC U4 

C A NEW ADJ CLASS 

C... ...RECORD A NEW ADJACENT FIELD NUMBER........ 

506 GO T0t517,5l7,5l7,5l7,5l7),NA0J 
NAD.)=NADJH 
GO TC518 

517 NADJ=6 > 

518 1FL=2 
MADJ (NADJ)=1B 
GO TC 114 

lltSET ADJCl.ASS 1FL=1 “ ' 

M. ^ 1. ■ • V ..fi? 
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C S £C n CN 2 

c *Sr4RT CaMPSLtNC TEHPQRARV STATISTICS * 

114 CD rO(5lO),ICRNUH 
11 ’J SUMeR=N 
NUriC-R=l • 

SUMER2=N*N 
ICRNUM=1 
JSL=J 
CO TC 102 

C... ...CCMTlNUe COHPIUNG TEMPORARY STATISTICS 

510 SU«L-R=SUMER+N 
NUH£R=NUPER*1 
SUMEH2=SUMER2+N*N 
GO TC 102 

C... ...IE GRAOIENr IS GREATER THAN OR EQliAl TO THE THRESHOLD 

C... ...CHECK IE PREVIOUS POINT A HAS A FIELD 

111 GO TC(5U,514J, ICRNUM 
GO TO IC2 

C YES .PREVIOUS POINT HAS A FIELD 

511 K6=0 
L0H=2S6 
LfiEAL=2!i6 

GO TCI221,22l,22l,221,13l),NADJ 
CO rC 22J 

221 GO rc(222,222) .NUHER 
GO rC 131 

222 LRCAL=5 
GO TC 144 

223 00 5112 1W=I>,NA0J ' 

3 IS L*M.U'J(lH) 

C... ...COMPARE LOWEST ADJACENT FIELD NUMBER 

',>04 IE (L-LOR)l 30,503,503 
150 KOaKH+l 
LQh = L 
LSIK81=10W 
LaLAILC^I 
GO TC504 

30J IFlLl'w.lt.LHCAULREAL-I.OH 
It) Hi 25 6 
■l.'J CON T! Nile 

L... ...tllH'AIl AOJACINT HELO (TRANSLATION) TABLE WITH IOHEST ......... 

l , ... A’.MACl'ir F !f I 0. ........ 

ill) 1 2H C'lil.KR 
K1--ISIK*)) 
l,’i> ( A(Kl 1 M (> AL 
nivTt''5U> 

I u u.n. 1 1 111 i 

1'.' liuai iCl* ■ 

0 ... HWIMIIM T li l.il NOMlIfH ro 2 'Hi 

1T<I> 111 14», 

MLAS;>J-1 

call fRANSI ITRRrR.MLAS.HHlN.MDUUn 

r,lj rC I4l 

144 CONfINUt ' 

l.A I LPLAL I..LKEAL 

512 Jl f)(J‘J-l 

C CODE FIELD SIUINCS IN THE DMA . 

DO 120 K'»=JSL,JLNO • 

12') Of:C(KO| = LMEAL 
,HADJ = 4 
IFL = 1 
tGRNLM=2 

NUPtLRFAL,U=NUH(LReAL, I )*NUMER 
SUM( LREAL, l )‘--SUH|LREAL. 1 l + SUMER 
SUH2 ILRfcAL.l ) = SIJM2ILRCAL,IM-SUMER2 
514 oeci J)*IDI hX(X,YI 
102 1FLG=IL’D 

C... ...HO PREVIOUS POINT A HAS NOT A FIELD 

C. A CLASS AT LINE END 

CO TC(51 1) , IGRNUM 
ID=3 . 

IC0=4 

CALL MASKOVIHZ, LIME, LIND 
CALL MASK0V(N7,0GC,DCCl) 

HRITE(PUIN,1010I (OECll I ), l»l,KDI 

101 continue 

HLASaJ-1 




C GO ro TftANSLATtON SUOROUTINE. 

CALL TRANS< lERR.K,MLASiMDIN|MOOUTJ 
lAi wruretLPfiooi) 

DO 130 1=6. ICO 
HUM=NUH( 1.1) 
lF(MO(>.EC.O)CO T0777 
AVE=SUMI I, n/HUM 

SO = tSUM2U ,n/NUKt I ,U-AVE*AVfi)**.5 

C... ...LIST STATISTICS Of- PERMANENT ACJACENT FIELDS.... 

no WRITEILPdOOZ) I.NUHt I . II . AVE. SO. SUHl 1 . 1 ) . SUM2 II . 1 ) 

CALL CLSAPE 
C FORMATS 

2000 FORMATC THERE ARE '. 15, • ELEMENTS IN CLASS*,t5,10X.I5,*A0J CLASS*) 

2001 FQRMATI' DEC',3UA) 

1007 FORMATI'O',' MAIN PGM DEC*/(32IA)) 

1015 FORMATI'O*, •IFLG='(13,' GRAD 0IR=*,I3| 
lOlO FORHArl20UOOA21 ) 

1000 F0RMATI6I5) 

1001 FURMATI • I' , lOX, '♦ AREA ♦ NOM OF PIXEL ♦ AVERAGE ♦ STANDARD 
1 DEV ♦ SUM • SUMSQR **) 

1002 FORMATI'O*. 10X,'7 '.lA,* ♦ '.IIZ,' * '.Fa.A,* • •.F13.A,* 

C* ' , III, ' ♦ • , I lA, ' ♦' ) 

i:'-: roR.-.A7!i;oi2coAi)) 

2006 FORMATC START CnilNTING A CLASS A* POINT', IS,' AND LIME', 15) 

20U4 FORMATC OH ODV CHOOY OHDOY QHDOY OH BOY', 15) 

2007 FORMATC NAO J = ' , 1 5 . lOX , • AD J CLASSsIB^* , I 51 

2008 FORMATC FINISH A CLASS', 2110) 

777 STOP 

END 

C 

C FUNCTION FOR COMPUTATION OF DIRECTION OF GRADIENT VECTOR 

FUNCTION lOIRXIX.Y) 

IF( ,NCT. IX.EQ.O.ANO.Y.eQ.O) )G0 TO 5 

IDIRX=0 

RETURN 

5 IF(X.EC.Q)X=>.001 
RAT=Y/X 

IF( .NOT. IRAT.GT.. A1A2.AND.RAT.lt. 2. AIA2) IGO TO I 

IDIRX=2 

RETURN 

I IFI.NCr. (RAT. lt. -.A1 A2. AND. RAT. GT.-2.AlA2nG0 TO 2 
tOIRX=A , - 
RETURN 

Z KAT'AQSIRAT) 

ID1RX=3 

IF(RAT.LT.2.A1A2)G0 TO 3 
RETURN 

3 tOIRXsl . 

RETURN 

END 

A** *• A* *t *>***« ******* ************************** **#**«•*«*•** 

c • 

C rKAASLAriCN ROUTINE 

SUHROUTIME TR ANSI I ERR , K.MLAS ,M0 1 N.MOOUTI 
C me US THE LOWEST CLASS NOT ALREADY PERMANENTLY ASSIGNED 

!HreCER*A NUN,SUM,SUM2 
INrEGER*2 DECU102A),D6C,LA 

LUHMUN NIiMI 2S6,P),SUMI2S6,2I,SUM2(256,2) ,IBC,ICD,lt ,LA(20Aai 
COMMON OfcCll02A).LCI256| 
logical*! CI20A8) 

LOlHvALENCEICm ,DECI 1)) 

LUC.ICAL*! 1C.A(S12)’,L0(512) 
lMrU‘.lK *2 0 ( 266 ) 

KOUIVALIMCE (OFCll 1 ) ,LCI II I 
ECUIVAISNCE (LAlK.AIll) 
tUUlVALlNCfc noil). Dll)) 
lU M =1 : 

01 .’ I =2 
HI n = T 

IX A) A • . 

(X ’.1 ■ • 

HI Hi NO MOn 
It I •• 

C 'A 


origin 
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wi' ( 1 1 (r t'. icB, n ,MLAS 
N(.li I . 

li'0& I lIKUrtl (•()•,’ ICI)-' me ELtfUNI 

WKftI (L»‘.100«HI,A( II, ICHI 

ICO/. IfJltf/AI ('(>*, I fAliLC LAin BCPOrtt JRAN5LAri0H‘/l*0«, 1214)1 

iarrsr«ii!C-i 

tm 2t>6 I so, ICO 
NUM( 1,71=0 
5UM|I,7|sO 
706 SUM? (1,7 1=0 
lMsICD-1 

C KtmiCi. CVEKALL ADJACENCY TABLE TO LOtaEST ADJACENT FIELDS 

DO 700 1=6, IFX 
H= I 

301 L=LA(M) 

IFIL.fcQ.flGO TC 300 

M=LA(LI 

00 TO 301 

300 IF||.-NCRT)25I,25l,2S2 

251 N0G-D(LA(L)» 

GO TC 240 

252 inO=IHCtl 
NGRT=L 
«UG=IBG 

C KCDUCE TEMPCBAftY STATISTICS INTO 

C permanent ADJACENT FIELD STA T I ST ICS . . . . 

240 NUHIMBG,2I=NUMIHBG,2I*NUM( I, 1) 

SUM(.‘'BG,2» = SUH(KBG,2)+SUM( I. n 
SUM2{P1)G,2I=.SUM2(HBG,2)*-SUM2I I, II 
NUMtI,ll=0 
SUM! 1,11=0 
SUM2II,ll=0 
D(n=MBG 
250 CONTINUE 

C.3«3«««* «««««« **«««;*•* A* •«****•*««»• 

C I6G IS NOW THE NEW LAST NUMBER FOR CLASES OR NEW (ICO) 

C LA IS NOW A VECTOR OF TRANSLATIONS FOR CLASS ADJACENCY 

C COMBINE OLD H 

C CCMBINE OLD NUM INTO NEW NUM.ISUM, SUM/) 

C*««444^A*4«*****4«*****A******44**««*«4*4*44****«**«««******«*«*«**«« 

IENO=II-l 
UU 257 1=5, iUG 
NUM( I, 1) =NUM(1 ,21 
SUM! I, IMSUMI 1 ,2) 

257 SUM2 (1 , 1 )=SUM2( 1,2) • 

CALL MASMQV(512,L0,LA) 

DO 254 1=2,512,2 

254 LO(I/2)=LO( II' 

WRITEILP*1005)(L0( I),I = 5, ICO) 

IC05 FORMATI’O', 'TABLE LO(IJ»/l'Q«, 3214)1 

c*.!; ."POT,*.".; ; , Trn,!.n, IC,-:,!) 

c 

c 

ICD=IDGtl 
HRITE(LP,lOll)KD 
toil FORMAT! 'O' , 'KD= *,I3) 

K0=K-2 

KE=2*K0 ' ■ 


HRITEILPflOOO) 

DO 260 KI=l, lENO 

C... ...READ temporary FIELD NUMBERS FROM DISC... 

REAOIMOIN, lOlOMCECKI) , I = l,KO) 

DO 250 M=2,KE,2 • 

250 LCIM/2)=LCIM) 

C TRANSLATE TEMPORARY FIELD NUMBERS TO PERMANENT FIELDS.......... 

CALL TRNTAB(0,KD,DECl,I6R,l) 

C... ...WRITE PERMANENT FIELDS ON TAPE 

260 WR[TE(HDOUT,1003MLCn),I=l,KD) 

M=M01N 
M01N=M00UT 
MOOUr=M , 

icoa FORMATCO' ,3214) 

10G-) FORMAT (• ','TRANSLATE.D DECl') 

H=2*--MLAS 
00 150 1=2, M, 2 
150 1,0(1/21=0(2) 


CALL TRNTAB(0,MLAS,0,tER,l) 
DO 151 1=2,M,2 
151 LC( I )=LUU72I 

iF( rco.EC.256) reRR»i 
ICOJ- F0R*/AT(1C:0(200A1)> 

ICIO FORHAT(20t 1C0A2) ) 

101‘) FORMAT t' ',32141 
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RETURN 

END 
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t-TEST COMPUTER PROGRAM 


CO>«PILER OPTIONS - NAPE- PAINtCPT>COi LiNECNT«6C( SlilE-OCOOK* 

SOURCE , BCOtNCL I S rtNOOECK.LOAO»NOPAP,NOfcD I T f NOlOtNOXREP 
DOUBLE PRECISION OEGF*. NCLEt SINA,COSA,FACTl,FACT2. JFACT.KPACT i 
♦ANCLEP,P1,TVAL»DVALUE»S1GLEV 
INTEGER F(2) .CI2) 

OIPENSION FMT(20I» FMHAOi. FM0I40) 

INTEGER STCLAS»OUTPT,STCLA 
DATA ICR.LP, 1. J/05*06,l,2/ 

DOUBLE PRECISION NUK( 2 ) • SUHXI 2 1 « SUHX2I 2 ) t AVE 1 2 1 f PVAL 

PI-3.iAlS926S358979328'<62A33B 
60 FCRPAT{4I5i2F5.2l 

62 F0RMAT(2CA<i) 

63 F0RPATI'2*> 

b<* FORPATI • ,100X,F12.6,F5.1I 

READ(ICRi60)NFEAT»NCLASS. INPUTvSTCLAStSIGLEVtOVALUE 
READI ICR.62)FMT.FMI,FH0 
00 S2 II-1*NFEAT 
STCLA«STCLAS 
V.RITE(LPt63) 

I“l 
J*2 

READ! INPUT .F HT ) I Ft K) tCI K I .NUH( K ) « SUHXI K I »SUHX2( K I »AVE(K ) .K«1.2I 
DO 52 JJ*2.NCLASS 

CALCULATE THE DEGREES OF FREEOQH AND X VALUE FOR THE 

t-oistribution . 

CEGF=NUM(I J+NUM( JI-2.0 

TVAL=(SUKX( I l/NUMt n-SUHX( J)/NUHIJ)-CVALUE)/(IDSCRT(1.0/NUK( D^l.O 
S/NUHI J) ) )*DSCRTnSUHX2l I) 4SUMX2 ( J )- 1 SUHX t I )«SUMXII ) ) /NUM i I I - 1 SUHX ( 
iJ)tSUNXUn/NUM( jn/CEGFH 


IDaOEGF 

LB>H00(IUt2) 

IF(LB.EQ.l) GO TO 20 

C. EVEN CALCULATION. 

IF(ID.GT.IOO) I0«100 
GO TO 24 

C ODD CALCULATION. 

20 IF(ID.GT.99I 10*99 
24 OEGF»IO 

ANGLE*DATAN(TVAL/IDSCRT(CECFn ) 


GO TO 1 421. ID 

SINA«OSIM ANGLEI 

COSA*OCOStANGLEI 

lD»ID-2 

PVAL*l.O 

JFACT*l.O 

KFACT-l.O 

IF(ID.LT.2I GO TO 30 
DO 26 IFA*2, 10,2 
FA- 1 FA 

GO TO 1251, LB 
JFACT*JFACT*(FAtI.O) 

GO TO 29 

25 JFACT-JFACT*(FA+k.O) 

29 KFACT*»(FACT*FA 
ANGLEP-l.O 
DO 27 IRP-lr IFA 
27 ANGLEP*ANGLEP*CQSA 
FACTl-JFACT 

ORIGINAi; PAGE IS 
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FACT2«KFACT 

CALCULATE THE SUNIFICANT LEVEL OF THIS SOLUTION. 

GC TC 132), LB 

FVAL-PVAL* ( ) r-ACTl/FACT2) •ANCLEP) 

GO TO 26 

32 PVAL-PVALftlFACT2/FACri)*ANCLEP) 

26 CONTINUE 
30 PVAL-PVALPSINA 
GO TO (36), LB 
33 GO TO 64 

34 PVAL>PVAL«fl2.0/Pl)*C0SA)tn2.0»ANCLEI/PI) 

GO TO 44 

42 PVAL«2.0*ANGLE/PI 
44 IFCPVAL.CE.SIGLEV) core 48 

C SIGNIFICANT LEVEL TO SMALL SUM TOGETHER CLASSES ANO RETRY. 

NUM( n-NUPintNUMIJ) 

SUNXm-SUHXID^SUPXIJ) 

SUHX2I I )«SUMX2I 1 )«SUMX2IJ) 

^R1TE(LP,FPI )F(I ),C( D.CIJI.AVEI 1),PVAL 
URITE<LP,64)TVAL,0ECF 

46 REA0llNPUT,FPr,EN0>S2)F( J ) , C| J ) , NUH| J ) , SUNK U I » SUMX2 1 J ) ,AVE IJ) 

CO TC 52 

46 LKlTb(LP,FNO)F( I l,CI I ) ,STCLA,AVEM ItPVAL 
mkITE<LP,6h) TVAL iOEGF 
STCLA»STCLA*l 
j;0 TO (49,50), 1 

49 I«2 

J«l . 

GO TO 46 I 

50 1*1 
J*2 

GO TO 46 
52 CONTINUE 
STOP 
END 





